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Chapter 1 
Introduction 



It goes without saying that an efficient method of performing accurate non- 
pcrturbative calculations in quantum field theory would be extremely useful. The 
method should be able to produce a sequence of approximations which can to be 
seen to converge. To be really useful, the approximation scheme should also be 
applicable when there is no identifiable small parameter to control the approxi- 
mation, and hence allow us to reach the areas of greatest interest in theoretical 
physics, eg the strong hadronic physics. In this thesis we outline a promising 
approximation scheme based upon the exact renormalization group. We start 
this introduction with a discussion of effective field theories. We then move on to 
discuss how we came to decide our method of approximation was a sensible one 
and discuss some of the problems posed by critical phenomena. 



1 



1.1 Effective Lagrangians 



The standard model of elementary particle physics has enjoyed considerable suc- 
cess. It can neatly account for most of the phenomena that are observed today to 
a high degree of accuracy. It is based upon the gauge group SU{3)c ® SU(2)l ® 
U{1)y- The 5'f/(3)c deals with the strong interactions that bind the hadrons to- 
gether , whilst SU{2)l ® U{1) deal with the electroweak sector. The Lagrangian 
density of the theory is, 

^SM = ^kin + ^yuk — ^Higgs + C-e (1-1) 

where Ckin represents the kinetic terms of the Lagrangian density and Cyuk are 
the Yukawa couplings which ultimately give the particles their masses. The 
term reflects the non-trivial vacuum topology of a four dimensional non-abelian 
gauge field theory. Perhaps the area of biggest speculation is that represented 
by Vniggs- It has long been postulated that particles gain their masses when the 
Higgs field , H, gains a vacuum expectation value, < H >= v. When this occurs 
the SU {?>) c® SU {2) L®U {1)y is spontaneously broken down to SU {?>)c ®U {1) em 
resulting in both fermionic masses (via the Yukawa couplings to the Higgs) and 
generating masses for the and Z*^ . This method of spontaneously breaking 
the symmetry has long been supported, although other competitors have existed 
(eg top quark condensates techni-colour, extended techni-colour [|, |] etc). 

Although the standard model has enjoyed considerable success there have been 
innumerable attempts to go beyond this theory, that is to extend it. One of the 
most promising early attempts was made by Georgi and Glashow who, based 
upon arguments of aesthetics and technical considerations, embedded the stan- 
dard model gauge group in a larger Lie group SU{b) and hence unified the three 
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gauge groups at a very high energy scale A, giving birth to the concept of grand 
unified theories. Since then considerable effort has been placed into investigat- 
ing such theories. For example, supersymmetry was introduced to cure what 
is known as the hierarchy problem (where the low energy SU{2)l x U{1)y 
breaking doublets receive radiative contributions to their mass of the order of 
the GUT scale). Other authors, interested in the lack of right-handed neutri- 
nos, have investigated left-right symmetric Pati-Salam models based upon the 
gauge group SU (4) ® S'f/ {2)i®SU (2)/j (and at the same time introduced a fourth 
'colour' !). Others interested in gravity have introduced super-gravity models and 
super-string models. 

The key point is as follows:- the extended theories always assume that the new 
theories only become important at high energies. To be more precise, if the 
extended theory is based upon a gauge group Q, at some energy scale A, then at 
some lower energy scales the gauge group is broken down to the standard model 
gauge group. Two important points now arise. Firstly the low energy physics 
depends on the high energy physics only in a limited way P, ^ |10|- That is, 
we are largely blind to any fundamental theory (ie a theory valid for all energy 
scales), if such a thing exists, and if it does we can only see small 'windows' of 
its effects on low energy phenomena. For example, it is usual for a grand unified 
theory (GUT) to predict proton decay. They also predict that the proton has 
a very long life time and is to all intents and purposes stable, eg the simplest 
SU{5) GUT predicts a lifetime of about 3 x 10^° years. Proton decay, despite 
being extremely rare, is one of the few windows open to us to look at the physics 
at the GUT scale. Of course this extremely long lifetime is important for our 
existence — if the proton decayed at short intervals then it would be unlikely that 
we would see a stable, evolving universe around us today. The weak effect of high 
energy physics on low energy phenomena is also important from a practical view 
point. If someone decides to investigate a new extension to the standard model. 
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then it must be checked whether the predictions he within the regions allowed by 
experimental results. For example, extended techni-colour was largely rejected 
by it predicting values for the so called S,T and U parameters that lie outside the 
bounds dictated by experiment. However, to compute results in these extensions 
requires the computation of a huge number of Feynman integrals, which leads 
to an arduous task. It is much simpler to realize that at the end of the day 
all that is required is to compute the values of a number of parameters in an 
effective theory, and thus to reproduce the result in a far more efficient manner. 
The point is that effective theories represent our best chance to systematically 
organize the experimental data and hence increase any chances of producing a 
worthwhile result. 

The second problem is a calculational one. We now have a problem with several 
mass scales:- in addition to the low-energy physics we now have a new large energy 
scale A. This vastly complicates the calculations. For example, within perturba- 
tion theory beyond tree level, internal propagators will have several mass scales 



involved and this leads to complicated Feynman parameter integrals ITT], |T2 



Perhaps a better way to proceed is to introduce the concept of an effective La- 
grangian. Instead of taking a theory to be valid over all energy scales, we instead 
assume that the theory is valid only over a limited range of energies, without 
requiring a detailed knowledge of the physics outside its range of validity. For 
example, the physics of soft pions scattering strongly below the chiral-symmetry 
breaking scale is described well by chiral Lagrangians and the standard model 
below the scale of grand unification. It appears that to some extent all the the- 
ories that we know are described by effective field theories up to the scale where 
some new physics occur. 
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We will no longer assume that the theory is valid for all energies and introduce 
an overall momentum cutoff Aq. Below Aq the physics will be described by a very 
general Lagrangian with an infinite set of couplings. Using 'naturalness' these 
couplings can be expected to be of order unity at the scale Aq. That is we will 
write, 

C,ff = Co + ^Ci + ^C2 + --- (1.2) 

where Co contains operators of canonical dimension < D, where D is the dimen- 
sion of space-time, and (for n > 0) contains operators of canonical dimension 
n + D. That is, Cq contains what are usually referred to as the renormalizable 
operators and the (for n > 0) contain the non-renormalizable operators. 

It may seem that we have lost all predictive power as our theory now contains an 



infinite number of non-renormalizable interactions, but this is not the case [121: 



if we know the underlying theory at high energy then we can calculate all 
the non-renormalizable interactions. 

Looking at expression ( |1.2D we see that it looks like the effects of the non- 
renormalizable interactions are heavily suppressed, unless the interactions 
lead to a divergence of a sufficiently high degree, which then overwhelms the 
suppression factor. For example, consider the six point coupling shown in 
figure (1.1a). If we consider the two loop graph in figure (1.1b) then we see 
that this will give a contribution of order Aq, which will cancel out the 1/Aq 
factor and lead to a term of order Aq. It again looks like we have lost all 
predictive power. However, the two loop graph will only alter the coefficient 
of the two point function (ie the mass term), which is a renormalizable 
term. Hence we can renormalize the two point function to remove this 



extra divergence. This argument can be shown to generalize [13 
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(b) 



Figure 1.1: Six point interactions in a scalar field theory. 



The fact that we can still calculate in an effective theory can be seen in one of 
the simplest of theories, quantum electro-dynamics. QED has managed to make 
spectacularly accurate predictions at energy scales from eV up to a few MeV. For 
example, the anomalous magnetic moment has been calculated to several loops 
in the effective field theory of QED. We say effective because corrections from 
the 'full' theory will be more important at higher energies, including the effects of 
QCD and electro- weak theory through loop contributions, plus contributions from 
any theory defined at a higher energy scale still. The suppression of these higher 
energy scale theories still allows QED to make accurate low energy predictions. 

We will look at the second point above more closely. We will define what is 



known as effective Lagrangian flow M, IT3|, IT5[. If we wish to look at the physics 



at some energy scale E <^ Kq, instead of using the full 'bare' Lagrangian with 
cut off Aq, we could lower the cutoff to a lower scale A ~ i?. To do this we will 
need to allow the couplings to flow so that the physics (ie the partition function) 
remains fixed for all the low energy processes. We will show later that the action, 
Sa = J d^xC, will evolve according to an equation of the general form (see |1^ 
and references therein), 

A^-ISl (1.3) 
6 



where Sa is known as the Wilsonian effective action. 



We no longer need to use the full theory at these low energy scales. Instead we 
use the Lagrangian defined at A to calculate physical quantities at the scale E. 
We define the effective theory to be renormalizable if we can calculate physical 
processes up to errors of order E/Aq, once we have determined a finite number of 



coupling constants at some scale A E [1^, |T5|, |T3- We call these couplings the 



relevant couplings and the other couplings the irrelevant couplings. We can now 
describe the low energy theory to a given order of accuracy without any recourse 
to the full theory or an infinite number of parameters. The flow of the couplings 
in the perturbative case is shown in figure (|1.2| ). Each trajectory or coupling flow 
corresponds to a particular choice of bare couplings at scale Aq. As we evolve 
the couplings down to the scale A we flow to a sub-manifold of the coupling 
constant space with a dimension equal to the number of relevant operators. This 
manifold has a thickness of order A/Aq (as we may expect considering that the 
effective theory is renormalizable, as defined above). We see that once the relevant 
couplings are known we automatically know the irrelevant ones to an accuracy 
of order A/Aq at the scale A, cf figure (pT^). 



Looking at figure ( |1.2| ) more closely we see that at low energy scales we have 
what is known as a 'self-similar' evolution — once we know the Lagrangian at 
low energy scale A, then the Lagrangian at a lower scale A' will look very much like 
the one defined at A, except for a change in the values of the relevant couplings. 
Now suppose we let the overall cutoff Aq tend to infinity, Aq ^ oo. We will no 
longer have a self-similar evolution, but will now have an 'exactly self-similar' 
evolution. That is lowering the cutoff will no longer affect the values of the 
irrelevant couplings (expressed in terms of the relevant couplings), and that the 
flow of the irrelevant couplings has reached a flxed point. This concept is known 
as the concept of universality. The analysis of flxed points and their classiflcation 
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Figure 1.2: The flow of the irrelevant couplings against the relevant couplings as 
the cutoff is lowered from Aq to A. The couplings flow to a manifold of thickness 
A/Ao for A < Aq. 



is central to the understanding of what possible theories could be relevant to 
the description of phenomena at present energies and is hence central to the 
understanding of what phenomena can be understood in terms of quantum fleld 
theory. Indeed, as Weinberg conjectured, any quantum theory that is Lorentz 
covariant, unitary and satisfles cluster decomposition is bound to look like a 
quantum field theory at sufficiently low energies |TB|, O . The illustration that 



this is true was based upon an extensive analysis involving the renormalization 
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group and fixed point analysis |]T5 . 

This is a convenient point to mention that so far all four dimensional field the- 



ories have been based around the Gaussian fixed point This consists only 
of the part quadratic in the fields and as such describes free, non-interacting 
field theories. As it is free it has no anomalous dimension. It is around this 
point that perturbative proofs of renormalizabilty and perturbation theory are 
based. It should be realized that we say based. The actual fixed point describes 
a non-interacting theory, whereas if we stay close to the fixed point, as in per- 
turbation theory, we can have a theory defined in perturbation theory with an 
interaction. We can perform the analysis of whether an operator is relevant or ir- 
relevant according to the dimension of the operator. Operators of dimension less 
than D (the dimension of space-time) are referred to as relevant, those of dimen- 
sion greater than D are known as irrelevant. However, it is usual to classify the 
terms in a Lagrangian due to their canonical dimension: we say that operators 
of canonical dimension less than D are relevant and those of canonical dimension 
greater than D are irrelevant. Of course it is well known that the dimension of 
true importance is the canonical dimension plus an 'anomalous' dimension due to 
quantum corrections. Hence, it may be possible for relevant operators, according 
to their canonical dimension, to become irrelevant, and vice- versa. Assuming that 
the anomalous dimension is small, we see that relevant operators correspond to 
power counting renormalizable interactions and irrelevant terms to power count- 
ing non-renormalizable terms. Therefore perturbative renormalizabilty amounts 
to the statement that quantum corrections grow at most logarithmically, so the 
anomalous dimensions remain small [13, 15 1. This is the case when we consider 
the Gaussian fixed point, provided we don't tune the irrelevant couplings to un- 
naturally large values. This is the theoretical basis upon which power counting 
renormalizabilty is based. 
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Some operators have zero canonical dimension, eg / d^x in D = 4. Oper- 
ators with zero dimension are known as marginal and perturbation theory with 
these interactions hence become important in their classification. That is, we 
can use perturbation theory with these operators to calculate (3 function of these 
operators. For example, it can be shown that / d^x (f){xY in D = 4 is in fact an 
irrelevant interaction [^, so it is not possible to define an interacting scalar 0^ 
in four dimensions with Aq = oo. (In fact Aq must remain less than some finite 
value.) This is a tremendously important point as 0^ theories play an important 
role in the Higgs sector of the standard model and if we have no interaction how 
can we possibly hope to break the SU (2)^ ® U{1)y symmetry of the electro- weak 
sector? Searches for non-Gaussian fixed points have been taking place for a long 
time (see for an early review), but have so far proved futileQ. Of course if the 
cutoff is kept finite then a quick glance at figure ( |1.2|) will show that it is possible 
to have an interacting effective 0^ theory in four dimensions. The key point is 
that the irrelevant couplings are suppressed by a function of Aq, which vanishes 
as Aq — > oo, but can play a role if the overall cutoff is finite. 

We have now come to a point where we should make the above ideas more precise. 
We begin with a discussion of fixed points, before moving onto the more detailed 
areas of deriving and approximating renormalization group equations. We will 
find it convenient to use critical phenomena as an example in later chapters, so 
we end this introduction with a brief outline of the theory behind problems in 
critical phenomena. 



^Seiberg and Witten have recently managed to find a non-trivial fixed points in A'^ = 2 
supersymmetric theories pl| 
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1.2 Fixed points 



One may ask what happens to the action as we continuously lower the cutoff. The 
simplest possibility is that it ffows into a fixed point of the transformation. That 
is an action, S*, such that continued lowering of the cutoff leaves it unchanged. 
Such an action will satisfy, 

t[S*] = (1.4) 



Other possibilities exist, such as limit cycles and turbulent behaviour, although 
they are less interesting from a physical point of view (see [|I^ and references 
therein). 

One may further enquire about the stability of these fixed points. To do this we 
linearize about the fixed point by writing. 



S[^] = S*[^] + 6Sa[(P] (1.5) 

This then yields, 

A J^<^^a[0] = -Li6S^) (1.6) 

where L is a linear operator acting upon ^^a. We will assume that L will have a 
discrete set of eigenvalues Aj corresponding to a set of eigenoperators Oi. In this 
case can be expanded as a series in the integrated Oj's 

5S^ = J2h,{A)0, (1.7) 

Equation (|1.6|) then yields 

A^/i.(A) = -X^MA) (1.8) 
Solution of equation ( |1.6|) will yield the values of the eigenvalues A,. 
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We define the critical surface as the surface, at A = Aq, containing all points 
which will eventually flow into a fixed point. We can classify the operators ac- 
cording to their eigenvalue:- 

If A > we say the operator is relevant. The operator will flow out of the fixed 
point as A is lowered. 

If A < we say the operator is irrelevant. Successive transformations on these 
operators force the operator into the fixed fixed point. 

If A = we say the operator is marginal. We can't decide what happens to the 
operator in the linearized approximation and have to go the quadratic approxi- 



mation to gain further information pO . 



It is possible for a renormalization group equation to have several fixed points. 
In this case we define the stability of a fixed point by the number of relevant 
eigenvalues it has. 

Given the above definitions we see that the critical surface is locally spanned by 
the set of irrelevant operators. 



This analysis is more extensively reviewed in p2| and [p3 



1.3 Renormalization Group Equations and Ap- 
proximations 

So far we have discussed renormalization group equations only in a general sense. 
Eventually we will wish to build an approximation scheme based upon the renor- 
malization group. In this section we will make the ideas of the first section more 
concrete by deriving an equation for a scalar field theory. We will then proceed 
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to discuss how to approximate these equations in a sensible efficient manner. 



1.3.1 Why use the RG? 

Perhaps the ffist question to ask is why should we use the renormalization group 
at all? For example why not use a scheme based upon improving the ladder 
ansatz in Dyson-Schwinger equations? The answer is quite simple - we wish 



to preserve renormalizabilty. It was pointed out in that Dyson-Schwinger 
equations quickly run into problems with renormalizabilty once we go beyond 
the ladder ansatz. This illustrates an important point — in any scheme involving 
truncations perturbative renormalizabilty is not guaranteed. 

We see that renormalizabilty is an important attribute to consider in any approx- 
imation scheme and must be preserved by the scheme. Polchinski [|13| applied the 
renormalization group to provide a conceptually elegant proof of the perturbative 
renormalizabilty of 0^ theory in four dimensions. It is only necessary to show 
that the irrelevant operators at zero coupling remain so at small coupling, but 
this must be true because the right hand side of the flow equation is a smooth 
function of the coupling. Hence, it becomes clear, both intuitively and in detail, 
that truncations of the flow equations are perturbatively renormalizable. The 
renormalization group seems a safe place to start the search for an approxima- 
tion scheme. 



1.3.2 The Legendre Effective Action 

We could derive an equation for the Wilsonian effective action S'a[0]. However it 
makes more sense to derive one for the Legendre effective action with an infra-red 



cutoff A [1161, E4|, E5i. There are two main reasons for this:- 
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1. It will be shown that the flow equations for the (Wilsonian) effective action 
Sa have a tree like structure with one particle irreducible bits linked by 



full infra-red cutoff propagators [|T3], |T6|, |2^. It will also be shown that 



we will wish to use a momentum expansion in our approximation scheme. 
This tree like structure must be preserved by any momentum expansion. A 
momentum expansion corresponds to expanding the vertices of the effective 
action in the scale of external momenta, regarding this as small compared 
to A. In the sharp cutoff limit this will cause all tree like terms with internal 
propagators to vanish, as, by momentum conservation, the momenta flowing 
through this internal propagator will be of the same scale as that of the 
external momentum. Noting that this internal propagator is also furnished 
with a sharp infra-red cutoff A in the sharp cutoff limit, we see that the tree 
structure will be destroyed |jl6|, ^ . This would be too much of a mutilation 



of the theory as all tree level corrections to the theory would be discarded 
as well as any loop diagram with more than one vertex. If instead we apply 
a momentum expansion to the one particle irreducible parts of S\ then this 
problem is avoided. 

2. It has long been known that we need to preserve a field re-scaling invariance 



if we wish to calculate the anomalous dimension rj ^7\, |3^ . ie 

(1.9) 

should be preserved as an invariance of the approximation scheme. 

If we use a momentum expansion of the Legendre flow equations then we 
know that such a re-parameterization invariance is preserved with certain 
choices of cutoff function, whereas there is no known one for the Wilsonian 



effective action Sa, when a momentum expansion is used pT| 
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The theory relating these two different types of action has been extensively de- 
veloped in [ |T6[| . 



1.3.3 Deriving an RG equation 



Our starting point will be the following effective action 




We assume the above has some internal symmetry, eg an 0{N) symmetry. [No- 
tation. It is perhaps a convenient point to explain the notation that will be 
used. We will used a condensed notation whenever possible. Hence J0 = 
J2a=i I xcf)"" (x) (x) , where a is any internal symmetry index. Similarly prop- 
agators are regarded as matrices in position/momentum space and any internal 
index space.] We assume the above is regulated by an overall UV cutoff Aq, and 
that is a smooth infra-red regulating cutoff function of width e. We assume 
it has the following property, C~^(g^) = {q-^^j^ — 1)?^. is a smooth regular- 
ization of the Heaviside 6 function, of width 2e satisfying < ^^(g. A) < 1 for all 
positive A and q, and O^iq, A) ^ 6{q — A) as e —>■ 0. 

If we differentiate with respect to A we get 



dA 




OA 



0exp-(-0.C-i .0 + J.0 + 5aJ0]) 




(1.11) 



Expanding the above then gives 



dWf,[J] 
dA 



2\ 6J ' dA ' 6J 



1 (6Wa[J] dC-^ 6Wa[J] 



+ tr[ 



dC-^6^WA[J] 
~dA 6J6J 



) 



(1.12) 



15 



Again in the above we have suppressed all internal indices. The first term in the 
above represents the tree like terms to which we objected before. The equations 
are best appreciated graphically as shown in figure ( |1.3|) . This shows how the 
n-point functions Wn{pi, • • • , Pn) of Wa[<^] evolve. 




Figure 1.3: The flow equations for the vertices of the generating functional W\. 
The vertices are shown as open circles. In a sharp cutoff the black dots will 
represent a delta function restricting momenta to g = A. The sum occurs over 
all possible disjoint partitions of momenta. 

We can now transform this to an equation for Fa, the Legendre effective ac- 
tion, with an infra-red cutoff A. We define Fa, the generator of the one particle 
irreducible parts of the vertices, by 



-WA[J] + Ji 



1.13) 



where (j)'^ = ^jj- is the classical field. Then using (derived in the standard way 
from ([rT3| )) 

s^w, fc-^^^y ,1,4) 



SJ5J 

we can re-write the equation as 



OA 



1 

-tr 
2 



dC 



-1 



5^Fa 



-1 



;i.i5) 



Again all internal symmetry indices have been suppressed. This will be our 
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starting point for all the work from now on 



Note that as e ^ the ^%]r- will restrict momentum to the shell q = A. Therefore, 



any other C ^{q) will contain 6'(0)'s, which are ambiguous |jT6|, Extreme care 
must be taken with this limit. 



1.4 Sharp vs. Smooth 

It has come to the point where we have to decide what type of cutoff to use. In 
this section we will briefly describe the pros and cons of sharp pO|, ^ and 
smooth [m, M cutoffs. 



1.4.1 Sharp Cutoffs 

The first problem that arises with using a sharp cutoff is how to take the sharp 
cutoff limit as e — > 0. We mentioned above that any other term containing a 
C^^{q) becomes ambiguous in this limit. To circumvent this problem we separate 
the problematic terms by writing 



{p,p') = 7(p, A)(27r)^5(p + p') + mip,p', A) 



;i.l6) 



so that f [0] = 0, and drop the field independent vacuum energy term. We can 
then write the subtracted equation as 



dA 



1 

--tr 
2 



dC- 
~dA 



C 



-1 



7 



'.f. (l + {C-i+7}.f' 



;i.i7) 



The sharp cutoff limit can now be taken. This is not so straightforward as it 
seems as we need to be able to deal with the ^(0)'s that occur when we do take 
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this limit. The answer turns out to be by no means as simple as using the usual 
physics convention of setting 9{0) = |! To take the limit we use the following 
lemma :- 



Let f{6^,A) be any function whose dependence on the second argument remains 
continuous at A = p in the limit e ^ 0. Then 



\Mp,A)f{e,{p,A),A)^6{A-p) / dtf{t,p) 

<:-*0 Jo 



where 5e{p,A) = —-^9e{p, A). This is easily proved by noting that 



1.18) 



(5,(p,A)/(^,(p,A),A) 



dA Je,{p,A) 



dtf{t,A') 



;i.i9) 



A'=A 



and that the integral is a representation of a step function with height dtf{t, A') . 
(This now yields 6e{p, A)9^{p, A)"^ — A), which doesn't equal the — A) 

which some may have naively guessed.) 



Using (|1.18|) , we may now take the sharp cutoff limit in equation ( |1.1?| ) to yield, 



9rA__l f d^q 5{q-A) 
'dA ~~2J (27r)^ 7(g,A) 



f. (l + G.f) 



;i.2o) 



where we have defined, 



lim 

C-Hp,A)+7(p,A) 



G(p,A) 



9{p - A) 
7(P, A) 



;i.2i) 



which represents a sharply infra-red cutoff, full two point Greens function. 



Notice that this equation displays field re-parameterization invariance under F — 
A2f,7-.A27 and G^G/X\ 

We now have a sharply cut-off RG equation. However, we quickly realize that 
any attempt to solve any of the above equations by a direct numerical approach 
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will quickly grind to a halt due to the sheer complexity of the problem |T^. We 
are therefore forced to choose some approximation scheme. 



As a first attempt we could try expanding in powers of field and truncating at 
some maximum order ie we write (assuming a Z2 symmetry), 

M 

Y = Y^a4^' (1.22) 

This expression is then substituted into equation (|1.2CI| ) and the equation ex- 
panded up to a maximum power of 0. This results in relationships between the 
coefficients which can easily be solved. This approach has been extensively in- 
vestigated. However few people at first recognized certain problems with this 
scheme. 



Firstly as pointed out in |3g] the truncations at first seem to converge to an answer 
but stop converging after a certain value of M. The reason for this is that there 
are in fact only a finite number of true solutions 0], together with an 



infinite number of solutions with singular behaviour for some real value of 0. Very 
bad solutions will have singular field dependence close to the origin causing the 
coefficients of 0^™" in ( |1.22|) to diverge with m. Of course, the truncation for which 



the coefficient of the 2M + 2 vertex vanishes will therefore better approximate 
the the Taylor series of a non-singular solution. At first increasing the value of 
M will improve the approximation, by forcing the singularities further away form 
the origin. However, even the non-singular solutions have singularities, but for 
complex 0, at a radius of, say , |0| = r. Therefore the truncations can not be 
expected to converge to better results than would be obtained with 'moderately 
bad' solutions with singular behaviour at or beyond the value of r. Also spurious 
solutions are generated and no completely reliable method can be found to reject 
these 
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Taking the above into account we could try some sort of momentum expansion. 
It should be noted that (|1.20| ) contains 6 functions, which do not have an analytic 
expansion. For example, 

^(|p + q| - 1) = 0(p.q + p/2) = ^^(q.p) + ^^^^^'^Hq.p) (1-23) 

m=l 

where we have defined p = p/|p| and is the m*^ derivative of the 5- 

function with respect to its argument. The above expansion means we cannot 
expand in powers of but are forced to use a non-analytic expansion in powers 



of momentum scale, p = ^ypJ), instead |16, 26 . 



To lowest order in the approximation we discard all external momentum depen- 
dence and write 

rA = Jd''x!^^{d,<pf + V{<P)^ (1.24) 
Such an ansatz yields the following equation, 

^-(^)-"(-^) 

This forms what is known as the local potential approximation and has been 
extensively investigated by several authors |^ ^ |3^ . It produces good 

results, although it does fail to take account of any wave function renormalization 
and other momentum dependent effects. 

The problems arise when we try to go beyond this level of approximation. To 
get any further, we will eventally be forced to calculate certain averages, which 



appear only to be calculable, in a closed form, in certain truncations [^. As we 
already know that truncations do not work we are forced to look down another 
avenue. 
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1.4.2 Smooth Cutoffs 



We have rejected sharp cutoffs and we Icnow that truncations don't work. We 
also know that any attempt to use a momentum expansion (that is a derivative 
expansion) must preserve a re-parameterization invariance. We also know that 
we need to keep all powers of the fields involved as we cannot use truncations in 
the fields. To satisfy these requirements we will employ a derivative expansion:- 

d^xV{cp) + ]^{d,<pfK{<P) + --- (1.26) 



This is substituted into ( |1.15|) and the right hand side expanded up to a maximum 
number of derivatives. 

The momentum expansion corresponds to a local derivative expansion in the 
effective Lagrangian, with radius of convergence p ~ e (this arises from the 
expansion of terms such 6'£(|q + p|, A) with q h). Since these expansions are 
substituted back into the flow equations and averaged over p ^ A we must have 
e > A for convergence. Obviously the cutoff must have A > £, as otherwise 
there would be no suppression of low momenta. Hence we see that e ^ A and 
typically the expansion will converge only slowly, if at all. To maximize the 
rate of convergence it clearly follows that we should choose the width to be as 
large as possible, as then we will have convergence as well as suppression of low 
momentum modes. We see that we are forced to use a cutoff that has a width of 



at least A 16, 30 



To preserve a re-parameterization invariance we will use a power law additive 
cutoff 

2k 

^(^'A) = |^ (1.27) 
for k a non-negative integer. The re-parameterization invariance manifests itself 
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in the form of a scaling symmetry with a set of (non-physical) scaling dimensions. 
That is if we choose 



[q^] = [d,] = l, [4>] = k + D/2, [rA] = 0, (1.28) 

as the non-physical scaling dimensions, then the scaling symmetry of equation 
( |1.15| ) becomes apparent. We already know that we need to choose a cutoff width 



of at least A and that convergence will be quicker the 'wider' the cutoff function 
is. As we have the following identity, 

2(fc+l) 

^ (1.30) 



j^2k ^ q2(k+i) 



we see that it is beneficial (ie results in the widest cutoff width) if we take k as 
small as possible. As we also need the cutoff to regularize the theory it can be 
shown that k > D/2 — 1. Hence we chose k to be the smallest integer greater 
than D/2 — 1. For D = 3 this means we take k = 1. 

The study of the above approximation scheme will form the major part of this 
thesis. We will apply it to some simple problems in critical phenomena and show 
that it does indeed work. We leave further development of the scheme to a later 
chapter and we will review some of the concepts in critical phenomena required 
for later chapters. 
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1.5 Critical Phenomena 



In this section we only be concerned with the problems posed by second order 
phase transitions. These are transitions where there is a continuous change in the 
properties of the system from one state to another, eg the continuous appearance 
of a magnetization as a ferromagnetic material is cooled. (As opposed to first 
order transitions where there is a discontinuous change, see for details). Such 
phase transitions are well known, even in everyday life. For example, consider 
water boiling in a kettle - the liquid is changing from a liquid to a vapour. At 
a certain critical temperature this phase transition is second order. In fact, if 
we look at this transition more closely we see the phenomena known as critical 
opalescence, where the fluid takes on a milky appearance at the transition point. 
This happens exactly at the transition point and is due to regions the size of 
microns fluctuating coherently on a large scale. This illustrates two difficulties 
with critical phenomena, 

• we need to consider all length scales. 

• we need to consider long ranges. 

This makes theoretical work extremely difficult; we find it nearly impossible to 
deal with problems involving just three degrees of freedom, let alone one involving 
maybe hundreds of thousands. 

There are other problems that we would like to understand in second order tran- 
sitions. For example, we would like to know why seemingly separate physical 
systems display very similar critical behaviour, eg uniaxial ferromagnets and flu- 
ids near their critical points display similar types of behaviour, even though they 
are totally different physical systems. This phenomenon is known as universality. 
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In the next section we describe some simple phenomenology of phase transitions 
and then briefly outline a theoretical background. There are many references to 



the subject but the books by Zinn- Justin p2[ and Amit p3[, and the article by 
Weinberg are particularly good introductions. 



1.5.1 Scaling 

For convenience we will mostly use the language of ferromagnetic systems. 

The basic quantity of interest in statistical mechanics is the free energy, F . This 
is defined by: 

F = -kT log Z (1.31) 

where Z is the partition function of the theory. For example, consider a classical 
spin model, so that the partition function takes the form 

Z = JU dS,p{S,) exp -pniS) (1.32) 

i 

where (3 = 1/kT, k is the Boltzmann constant, T is the temperature, p{Si) is a 
spin weighting factor describing the local microscopic properties of the system, 
and Ti is the Hamiltonian given by, 

J2V,,S,S,+Y.HiSi (1.33) 

ij i 

The variable Hi in the above represents an externally applied magnetic field. Such 
models have achieved great success in describing simple ferromagnetic systems. 
The statistical average of a quantity X is given by 

<X>=\ j X{dSip{Si)Xie^^-(37i{S) (1.34) 
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Using this definition we define the magnetization of a system to be the average 
of the spins 

M{H,T) = V-^^a^ < Sn> (1.35) 

n 

where a is the lattice spacing. In what is known as the thermodynamic hmit, 
where N ^ (yo, we can describe the magnetization as 

M{H,T) = -^F{H,T) (1.36) 

This can easily be seen by realizing that taking a derivative with respect to H 
brings down a Si into the sum. It can happen, that even in the presence of 
zero external magnetic field, below a certain temperature Tc there is a non-zero 
value for the magnetization. When this occurs there is said to be a spontaneous 
magnetization and the temperature at which this first occurs is called the critical 
temperature. 

There are other quantities of interest. The susceptibility is defined as 

X(H,T) = ^Wa'^< (Sr.- < S >)iSm- <S>)> (1.37) 

and this represents the response of a system to a small applied magnetic field. 
Also of interest is the specific heat which is defined as 

^--^df^ (1-38) 

Various physical quantities diverge as we approach a second order critical point. 
For example, the correlation function of two fluctuating fields behaves as 

9nm =< {Sn- < S >){Sm- < 5' >) >~ exp(-|a;n - Xm\/0^'' SiS X ^ OO 

(1.39) 
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We call ^ the correlation length. As we approach the critical temperature the 
correlation length diverges as follows 



e~ (1.40) 

{T.-T)-"' T<T, 

We call u a critical exponent. Other exponents can also be defined. For example, 
the susceptibility diverges as 

(T-T^-^ T>T^ 
X- , (1-41) 

(T.-T)-^ T<T, 

and the specific heat acts as 

C~ ^ ' (1.42) 

(Te-T)-"' T<T, 

Below the critical temperature a spontaneous magnetization appears which be- 
haves as 

M ~ (T; - Tf (1.43) 

If we are exactly at the critical point, T — Tc, then the correlation function will 
act like 

< {Sn- < S >){Sm- < S >) >r^ _ ^ |£)_2+r, a < |a;„ - Xm\ < oo 

(1.44) 

Notice that we have lost all functional dependence on a fundamental length scale 
— this reflects the fact that at criticality the system is scale invariant, and is also 
consistent with the appearance of long range order. The exponent rj is known as 
the anomalous dimension. 
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Various relationships exist between these exponents defined at the fixed pointQ, 
known as scahng relations. These mean that you only need to know two of the 
above critical exponents to calculate the rest. For example, 



2/3 + 7 = 2 -a (1.45) 
7 = u{2-ri) (1.46) 



are the two of the best known. 

These relations are easily reproduced. As an example consider the relation ( |1.46| ) 



To prove this, first consider the correlator of two spins (|1.44|). We see from 



equation ( |1.44|) , that when we are close to the transition point that the correlator 



2J_) 

has dimension —D + 2 — 77. Hence, if we consider ^y- 9ij then this will have 
scaling dimension 2 — 77. Therefore, as the correlation length sets the basic scale 
of the system at criticality, so that at criticality the system loses all dependence 
on all other length scales, we can determine the behaviour of the magnetization, 
specific heat etc, in terms of the behaviour of ^. This hypothesis is known as the 
scaling hypothesis . Using the scaling hypothesis, we see that 9ij satisfies, 

E 9^J ~ ~ - T,y(-'+^'> (1.47) 

near the critical point. However, a quick glance at equation (|1.37|) reveals that 
X = ^^J2ij9ij- Therefore, we must have 

X ~ (T - Tc)-^ ~ (T - T,y^-^+"^ (1.48) 
This then forces equation (|1.46|) to hold, ie 7 = z/(2 — t]). This is known as 



^The exponents discussed so far are defined by the behaviour at the fixed point. In addition 
to these there is also an infinite set of exponents characterizing corrections to the scahng laws, 
when the system is close to the fixed point 
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Fisher's scaling relation. It is also found that the indices are symmetrical about 
the transition point so a = a' , u = u' etc. The reader is referred to the extensive 
literature on the subject for further details. 



1.5.2 Critical Phenomena and the RG 

In this section we will indicate how to relate the above theory to the renormal- 
ization group. We know that near criticality the correlation length is very large 
and we need to consider a large number of degrees of freedom. Hence, if we wish 
to perform calculations near the phase transition and gain a theoretical under- 
standing of the problem, we need to find some method of reducing the number 
of degrees of freedom. The renormalization group does such a thing. If we con- 
tinually integrate out the fluctuations with the shortest wavelength we gradually 
reduce the degrees of freedom, making the problem more tractable. The key point 
relating this to the critical point is that at criticality all wavelengths are equally 
important. This is because the correlation length diverges and we therefore lose 
all dependence on any fundamental length scale (cf equation ( p..44| )). This also 



means that as we continually integrate degrees of freedom out of the theory we 
essentially see the same picture — ie we are at a fixed point of the transformation. 

To calculate the critical exponents we need to consider the behaviour near the 
fixed points. To aid this discussion we will re-write the equations in terms of 
dimensionless quantities. Suppose we are near a fixed point, but not on the 
critical surface. Equation (|1.8|) shows that 



~ A"^' (1.49) 
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As the hiS are now dimensionless we see that we have 



(1.50) 



Then, as the flow of the action wiU be dominated by the largest relevant eigen- 
value, we can write, 

SSA-^l^f^y^OL (1.51) 

where A^, is the largest eigenvalue and Cl is the operator associated with it. We 
see that the integrated operator Ol has dimension — A^, and is associated with a 
couphng fjL^^ of dimension A^. Now consider the situation at A = Aq. There will 
be an operator corresponding to the deviation from criticality. This will be the 
dominant operator in determining the statistical state of the system. Therefore, 
we associate this operator with O^. 

5Sao = (T - n)OL + 0{{T - T,f) (1.52) 

As we know that the coefficient oIOl has dimension Al, we see that the scahng 
dimension of the temperature difference is also A^. Therefore, as [^] = —1, we 
must, have 

^ - (T - (1.53) 

or, V = 

The value of v is therefore determined by the linearization procedure. The anoma- 
lous dimension is a property of the fixed point itself — in the full un-approximated 
equations the value of r] is determined by the field re-parameterization invariance. 
The scaling symmetry of the equations turns the renormalization group equations 
into a non-linear eigenvalue for ry, with only a few values of 77 leading to accept- 
able fixed point solutions. By preserving the re-parameterization invariance in 
an approximation scheme we can still determine 77 without recourse to any non- 
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physical arguments or parameters. We can now determine two critical exponents 
and therefore determine the others. Corrections to the scaling behaviour are 
given by the subleading eigenvalues. For example, the leading correction to scal- 
ing exponent, denoted hj u, is defined as being equal to minus the least negative 
eigenvalue. For example, close to the critical temperature we have the following 
scaling relation, 

^^\T-Tc\-'' + --- (1.54) 

However, there are also corrections to this expression, so further away from the 
critical temperature we have, 

^ ~ |T - T^]-" + a^\T - t/'^~^> + ... (1.55) 



1.5.3 Critical Phenomena and Field Theory 

We have now reached the stage where we should relate physical systems to field 
theoretical models. It should have been realized by now that there are two equiv- 
alent descriptions of critical phenomena. In the above we have sometimes used 
the language of classical spin like systems and sometimes we have found it more 
convenient to use a field-theoretical like language. There is in fact an intrinsic 
link between the two descriptions and we will make an attempt to briefiy describe 
the relation between the descriptions. Consider the classical spin model defined 
above by equation ( |1.33|) . If we assume that the lattice is hyper-cubic and that 
we only have nearest neighbour couplings then we can write 

J2 V^JS,S, =KJ2 (T.iSn+^ - Sn)' " 2dSl\ (1.56) 

ij n \ 11 J 

where /z runs over the directions of the lattice. In fact, if we consider a model 
where the spin is constrained by S"^ = ±1 (known as the Ising model), we can 
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describe this using a spin weighting function of (taking an appropriate choice of 
A and k), 

p{Si) oc exp -{kSI + XS^) (1.57) 
We can then write the Hamiltonian as 

n = Y: (K{f3) J2iSn+, - + f^{f3)Sl + X{(3)S'^ + H^S^ (1.58) 

n \ /J. ) 

[It is perhaps worth pointing out that other versions of the spin-weighting function 
( |1.571 ), with terms of higher powers in could also be used. However it turns out 
that such terms are irrelevant, in the sense defined above, and have no influence 
on the critical theory.] 

The above expression ( |1.58| ) for the Ising Hamiltonian should look familiar — it 



looks like a lattice regularized 0^ theory [45]. To see this take a massive 0^ in D 
dimensions, 

and place it on a hyper-cubic lattice by writing. 



X ^ Xn = an^e^ (1.60) 

(f){x) 0„ = (p{Xn) (1.61) 

D(p ^ l[d(j)r, (1.62) 

n 

df^Cp ^{(f){Xn + e^) - (f>{Xn)) = ^{(pn+fi- (Pn) {1.63) 

where a is the lattice spacing and is a set of D orthonormal vectors. Under 
these transformations the lattice regularized theory is 

S{<j), ^) = E j 4 E ^ ('/'n+M - + + ^K. + Jn<Pn] (1-64) 
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and the partition function becomes, 



Z{J) = ll f d(j)nexpS{(f),J) (1.65) 
where the functional measure D(f) is replaced by a finite dimensional integral 

Tin d(Pn- 



The correspondence between ( |1.59D and ( |1.58D should now be clear. So we see 



that there is a direct link between lattice theories and the Ising model. However, 
what we are really interested in is the critical theory. To clearly see the link we 
will have to show explicitly what happens at the phase transition and consider 
the ground states of the Ising model. 

Consider the potential of the model /i(/5)5'^ + A(/5)5'^, where for the moment 
we will ignore quantum corrections (ie we will use what is known as mean field 
theory [Q). For > we have a potential like figure (pT^). We see we have 
a ground state in which Sn = and 

< S >= ^Y.^n = for > (1.66) 

If ^{(3) < then we have a potential like figure (|1.4|b). Now the ground state 
corresponds to the Sn aligned with = w so that 



<^>=y^ forM/?)<0 (1.67) 

We see that the phase transition corresponds to the point where /i = 0. The 
temperature that corresponds to this is called the critical temperature. 

Now let us make the correspondence between the critical Ising model and and 
the critical lattice theory. To make direct comparisons we first move to non- 



32 



(a) /i > 



(b) ^ < 



Figure 1.4: The potential, fi{j3)S'^ + X{P)S'^, for the unbroken case yU > (a), and 
the broken case /i < 0, (b). 

dimensional variables. That is we write 

A = /a^-^ (f)' = ga^-'(f), and J' = ga^'/^+^J (1.68) 
the action then becomes 

S{<l>\ J') = -4 E U E(</>U. - <l^'rf + \^'<^'€ + + J'<^'n] (1-69) 

We see comparing equations ( |1.58| ) and ( |1.69| ) that we have the following corre- 
spondence, 

rn^a^ ~ (1.70) 

We know that as we approach the critical temperature /i(/5) — > 0. This forces us 
to fine tune a to zero as well, as the physical mass of the system m must be held 
fixed. Hence, we see that as we approach the fixed point we must also approach 
the continuum limit. 
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We should really have expected this on heuristic arguments:- we know that at 
the critical point the Ising model is scale invariant. Therefore to preserve this 
invariance we cannot have any massive parameter in the field theoretical model 
and so there cannot be any cutoff involved. We can also use this argument 
to show that the critical theory corresponds to a massless theory — any mass 
involvement in the critical theory would break scale invariance. We see that the 
critical theories correspond to massless, renormalized theories. 

We now have a definite link between field theories and the statistical mechanics 
of critical phenomena. There is one final missing link in our puzzle — how 
to relate the theory to physical systems. This is not as complex as may first 
be thought. For most simple systems it is possible to find local observables 
whose values depend upon the phase that they are in. For example, in the 
above the spin played the role of the order parameter, differentiating between the 
phase with a spontaneous magnetization and the one without. The relation of a 
physical theory and a theoretical analysis boils down to determining what are the 
order parameters and what are their internal symmetries. For example, the order 
parameter for the helium superfluid transition is complex having a continuous 
symmetry corresponding to multiplication by a phase. This can be shown to be 
equivalent to a theory of two real scalar fields. Therefore, we may use a field 
theory consisting of a real scalar field with a global 0(2) symmetry. 

We will primarily concentrate on systems described by 0{N) symmetric scalar 
field theories with the following action 

^ = /rf^a; 1^ (5^$^(x))' + ^$-(x)$-(x) + I ($-(x)$-(x))^| (1.71) 

As usual the index a refers to an internal symmetry under the 0{N) group. These 
models correspond to different physical systems according to the value of N (see 



1 2^ for details on how to derive these). 
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N=0 critical behaviour of polymers. This is only formally defined in the limit 
iV — > 0, as first noted by de Gennes p6| . 



N=l the liquid- vapour phase transition |47|, the alloy order-disorder transition, 



uniaxial ferromagnets. We see that here we have no internal symmetries 
involved, eg the liquid- vapour transition can be modelled by particles living 
on a lattice: allowing the occupation of each lattice site to be either be 



or 1 leads to a link with the Ising model ES 



N=2 He^ superfluid phase transition, planar ferromagnets. The first example 
was discussed above. The planar ferro-magnet will clearly have a symmetry 
in a plane. This then leads to the choice of a two component field with an 
0(2) symmetry between the two components. 

N=3 ferro magnetic phase transitions. Now we have a true three dimensional 
symmetry so a three field order parameter is chosen, with an 0(3) symmetry 
between the components. 

N=4 It has been postulated that this corresponds to the chiral phase transition 
for two flavours of quarks . 



The fact that one model for a particular value of can descibe several physical 
systems is the theoretical basis of universality. It should also be noted that it 
doesn't matter what the initial value of any coupling is, or the lattice spacing 
or even upon the shape of the lattice: once we have adjusted the temperature 
(or pressure, density etc,..) so that we are on the critical surface, we know that 
renormalization group transformations will take us into a fixed point, which will 
yield the same critical behaviour. In fact we can take this further and ask why do 
different materials have the same critical behaviour, eg in critical binary fluids, 
which correspond toA^ = l,D = 3,a mixture of aniline and cyclohexane shows 
the same critical behaviour as a mixture of triethlamine and water |5^. This is 
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because the critical behaviour is independent of the underlying structure of the 
physical system. At distances of order of the lattice spacing the detailed fine 
structure, such as the lattice spacing and the lattice structure, will be important. 
However when we use the renormalization group to decrease the amount of 'mag- 
nification' these details are washed out, as the details at the scale of the lattice 
will be averaged out [pll . 



1.5.4 Other methods of investigating critical behaviour. 

Before going any further we briefly consider other methods of investigating critical 
phenomena and calculating exponents. 

Original predictions of critical exponents were made using mean field theory in 
which a saddle point expansion is made about the classical minimum of the 
field equations. Exact predictions for the exponents are made. In fact for the 
0{N) models descibed above the predictions are independent of the value of 
A^. The main problems occur when we try to go beyond the tree level expansion. 
Below four dimensions infra-red divergences plague an expansion making it useless 

©El- 

To go beyond mean field theory Wilson and Fisher developed what is known as 
the e expansion. The dimension of space is taken as 4 — e and a double expansion 
is made in the coupling constant g of the theory and e. Thus, for example, a 
n-point function may be expanded as 

r"(pi,...,p„) = 5:r(:Ve^ (1.72) 

rs 

This coupled with a renormalization group analysis has produced some spectac- 
ular results. This method avoids the infra-red divergences by moving from a 
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gaussian fixed point to the non-trivial, infra-red stable Wilson-Fisher fixed point 
The e expansion has enjoyed considerable success and has made accurate 



predictions of the exponents. However it does have some drawbacks. For example 
the series it produces is only asymptotic and requires Borel resummation ||2^ . 



The N vector model has also been studied on lattices with nearest neighbour in- 
teractions. The results for the exponents come from analysis using high temper- 
ature series by different types of ratio method, Fade approximants or differential 
approximants. For large N , the expansion has also been used. However the 
large N models are unphysical and are only studied out of theoretical interest 
(the main one being that the N = oo model is exactly solvable). 

There are various other methods that have been used to calculate critical expo- 



nents which deserve a mention. The real space renormalization group method |52 
consists of considering spins on a lattice and performing blocking transformations 
to reduce the number of degrees of freedom, and then truncating the resulting 
expression at a certain number of operators. Similar in style are Monte-Carlo 



renormalization group calculations |5J, Here the path integral is ex- 



panded in a certain set of operators and a blocking transformation performed. 
A Monte-Carlo calculation is then performed to calculate what values the cou- 
plings associated with these operators should have afterwards. Finally there is 
perturbation in fixed dimension p2|. Developed by Parisi, this consists of calcu- 



lating the f3 function in perturbation theory at a fixed dimension. The resulting 
expression can then be re-summed to provide information about the zero's of the 
f3 function and hence information about the critical exponents. 
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Chapter 2 

The Results at Leading Order 



In the introduction we outlined why an approximation scheme based upon a 
derivative expansion of the renormahzation group equations is sensible approach 
to analytic non-perturbative calculations. In this chapter we will report the 
results of applying the derivative expansion at the leading order in the approx- 
imation to a theory with a global 0{N) symmetry, and calculate the critical 
exponents 77, u and cu. 



2.1 Deriving the Leading Order Equation 



Our starting point will be the flow equation for the Legendre effective action, 



OA 



-tr 



+ 



5(f)"-5(f)^ 



(2.1) 



where we have dropped the superscript c on the 0's, denoting the classical field, 
and re-introduced the internal symmetry indices. The trace represents a sum 
over the internal spin indices and an integration over momentum/position space. 
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As in [0 it is convenient to write the trace as an integral over momentum space 
and factor out the D-dimensional sohd angle: 



dA 



D-l 



dCiq,A) 



^C(g,A) dA 



5"^ + C. 



n -1 



(q,-q) (2.2) 



where VL = 2/[T{D/2){4:n)^/^] is the solid angle of a D — 1-dimensional sphere 
divided by {2tt)^, the brackets (■ ■ ■) represent an average over all directions of the 
momentum q, and the trace now represents only the trace over the spin indices, 
a and b. 

It is not yet clear how the anomalous dimension will be determined. At a fixed 
point we know the field scales anomalously as ~ A^0, where = \{D-2+'q) 
and 7] is the anomalous scaling dimension. Considering the defining expression 
for Fa, a[4>\ = — |0.C~^.0 — Wa[J] + J.0, dimensional analysis then shows that 
we require C to behave as follows, 



C{q,A)^ A'^-^Ciq^/A^) 



(2.3) 



for some C, if we are to have Fa independent of A as we approach the fixed 
point. From now on we will write C{q, A) as A''~^C(g^/A^) and drop the tilde on 
the scaled C . As we are only interested in the behaviour near or at fixed points 
this is a sensible definition and we can take to be a constant. This is another 
reason for choosing an additive cutoff. A multiplicative cutoff would not scale 



anomalously, as is required if we wish to reach non-Gaussian fixed points |30 



From this discussion we see that it will be particularly convenient to re-write the 
equations in terms of dimensionless variables. 



q 



0(Aq) 



Aq 



iq) 



(2.4) 



39 



We will also re- write the equations in terms of t = log(Ao/A) and rescale the 
fields and the effective action to absorb the factor of in ( p.2|) , 



r 



2C 



where we have dropped the subscript on F, and C is a normalization factor to be 
chosen for later convenience. Upon doing this we get, 



(2.5) 



-Ctr/ rfgg^-i 
Jo 



' q dCjq^ 
C{q') dq 



+ 2-7] 



6^T 



iq, -qj 



In the above A^p = is the field counting operator: it counts the number of 
occurrences of the field in a given vertex, and arises due to the scaling of the 
field in ( 12.4] ). Aq is the momentum counting operator plus the dimension of space 
D and arises through the rescaling of the momenta in equation ( |2.4|) . It can be 
represented as 

Operating on a given vertex it counts the total number of derivatives acting on 
the fields 0. 

Equation ( p.5| ) will be the starting point for all the work from now onwards. 
Notice that for the first time rj explicitly appears in the equation. 

We write r[0] as the space-time integral of an effective Lagrangian expanded in 
powers of derivatives, 

m = J d^x {v{<p\ t) + ^K{<p\ t) {d,rf + \z{(^\ t) {ro.rf + ■■] m 
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Each linearly independent (under integration by parts) combination of differenti- 
ated fields will carry its own general (t-dependent) coefficient. The global 0{N) 
symmetry forces us to choose the coefficient functions to be functions of 0^. We 
will require that the fixed point solutions for V,K,Z etc will be non-singular for 
all 0^, that perturbations about these solutions grow no faster than a power, and 
that K{0) ^ 0. 

The rescaling symmetry is made explicit by choosing the following (non-physical) 
scaling dimensions, as follows from ( p. 71) and the definition of F. 

[qf^] = 1 [^^'] = 1 [0"] = k + D/2 (2.8) 

[V] =D [K] = -2{k + 1) [Z] = -2{2k + 2) - D 

where k is the exponent in the cutoff function, C{q'^) = q^^. The expansion is 
performed by substituting equation ( p.7|) into equation (p.5|) and expanding the 



right hand side up to a maximum number of derivatives (in this chapter this 
will be no derivatives). The angular average in ( |2.5| ) can be easily computed by 
translating them into invariant tensors, eg (g^g*^) = q^5^^ /D etc. From now on 
we will specialize to D = 3. 

To lowest order in the expansion we drop all the derivatives from the right hand 
side of ( p.5| ). The coefficient functions K,Z, etc are then determined by linear 
equations given by the vanishing of the left hand side of ( p. 51) . Consider the 
equations for 

^K{<P\ t) + (1 + vWK'{<P\ t) + T^K{ct>\ t) = (2.9) 
We see that at fixed points, where d/dt = 0, this then predicts 

K{(j)^) oc (02)-r?/(i+^) (2.10) 
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where ' — To avoid K{(tP) being singular and ensure that i^(0) ^ we see 
we must have 77 = 0. Therefore K{4P') must be a constant. Using the scahng 
symmetry we can set this constant to be 1. 

If we now consider the equation for Z we get 

t) + 4>''z\4>\ t) + z{(t>\ t) = (2.11) 

Hence, we see that at fixed points that Z will satisfy, 

Z oc (02)-^ (2.12) 

Again to avoid non-singular behaviour we set the constant of proportionality to 
zero and get Z{(j)^) = 0. In general, considering the form of a general term in 
the expansion we see that a similar conclusion must hold for all other terms. A 
general term, will satisfy 

^//(0^ t) + (0'i/'(0', t) + n^H{(t)\ t)) + n9H{(t)\ t) - 3//(0^ i) = (2.13) 

where is the number of 0's occurring in the term multiplying if, divided by 
two, and rid is the number of derivatives occurring in the term multiplying H. 
Hence we see that at a fixed point, where rj must be zero, we have 

H oc (02^-K+«9-3) (2.14) 



We see that for terms of higher order in the expansion, where > 2 and ng > 2, 
that H will be singular unless the constant of proportionality is set to zero. 
Therefore, at leading order, the derivative expansion must reproduce the local 
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potential expansion at fixed points: 



a\2 



(2.15) 



The next step is to calculate the trace of the following expression, 



6"'' + C. 



iq, -qj 



5''^ + 2g2 + 20^^"' + gV'}] ^ (2.16) 



Working in momentum space and integrating by parts where necessary we have, 



5^m 



25"V + 40''0V' + q'^S"'' 



(2.17) 



The taking of the trace over the spin indices of ( |2.16| ) is easily accomplished 
by noting that an operator y45"* + B(j)"-(j)^ has eigenvalues A and A + with 
the eigenvalue A having multiplicity — 1, where N is the dimension of the 
internal symmetry space. Hence, provided that A ^ and A + (p'^B ^ 0, then 
+ has eigenvalues 1/A, with multiplicity A^-1, and 1/(^ + 0^5). 

Hence we must have tr[A5''* + = (A^ - 1)/A + 1 /{A + BcfP'). Therefore, 

taking the trace of equation ( |2.16| ) gives 



N -I 



+ 



1 



1 + 2V'q^ + g4 1 + {2V' + 402y")g2 + g4 



(2.18) 



The leading order equation has now become 



dv{<p\t) 



dt 



AC dqq' 



N -1 



+ 



1 



(2.19) 



l + 2V'(02,t)g2 + g4 1 + (2V'(02^t) +402-l/"(02,t))g2 + g4 



Performing the g-integrals, doing the trivial angular average and setting ^ to l/27r 
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yields the equation at leading order, 



dt 



^2 + 2V'(0^ t) + 402-i/"(02, t) 

(2.20) 



2 + 2V'{(P^,t) 



In the rest of this chapter we will report the results of the study of this equation. 



2.2 Fixed Point Solutions 



We are now ready to start the search for fixed point solutions of ( p. 20 ), that is 



solutions with d/dt = 0. At first sight it may seem that ( p.20|) has infinitely 
many solutions parameterized by the value of V"{0). In fact this is not the case, 
as only finitely many solutions do not end in a singularity |^0[ Of course 
this is sensible on physical grounds, as the fixed points correspond to massless 
continuum limits (ie second order phase transitions) with the prescribed field 
content. To see that we should only expect to see a discrete set of solutions we 
must consider the boundary conditions supplied to ( p.20| ). 

The only requirement we have imposed upon V{(j)'^) so far is that it is non- 
singular for all 0^. Considering the form of ( p.20| ), we see that this implies that 
either V{(f)'^) is trivial (the Gaussian fixed point), or that V(0^) must satisfy the 
following for large 0^, 



A(^^)'-.^(^S^^ + 0((,=)-) (2,21) 
for some constant Ay. 

If we consider the perturbations about this solution we see that at large 0^ the 
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perturbation must behave as a combination of exp(|[3OA„]^/^(0^)^/^) and (0^)^. 
The second perturbation merely alters the value of Ay, whilst we disallow the first 
as it is incompatible with ( p.21| ). We see that the solution space of fixed point 



solutions, defined for all (f) and satisfying (2.21), divides up into an isolated one 



parameter set, parameterized by the value of A^. 



The only other possibility that makes the two sides of (|2.2(]| ) balance is that V{(f) ) 



is only defined for (f) < (p^ and ends at a singularity, at (p = (p^ as follows, 

1 /27\2/3 . x4/3 

(0c - 0') (2.22) 



16/ 



where we have suppressed non-singular and lower order-singular parts. Notice 
that if either V or V where to diverge then it would not be possible for both 
sides of equation (|2.2CI|) to balance. We will disregard these non-physical singular 



solutions. 

The imposition of the 0{N) symmetry forces to be a function of 0^, so it is 
symmetric under the Z2 transformation 0" ^ —(f)"'. To ensure this, we need V{(f)'^) 
to exist at the origin and satisfy ( |2.20| ) there. Setting 0^ = in equation ( |2.20| ) 
gives us our final boundary condition, 

- 31^(0) = - . (2.23) 
^2 + 2V'{0) 

We now have a second order differential equation with two boundary conditions. 
Thus we expect at most a discrete set of acceptable solutions. In fact we only 
find two:- the Gaussian fixed point and an approximation to the Wilson-Fisher 



fixed point |[T4 |. 
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2.2.1 The Gaussian Fixed Point 



It can be quickly seen that imposing = and = that a trivial 
fixed point solution exists, 

V(^^) ^ ^ (2.24) 
This is the fixed point mentioned in the introduction. 



2.2.2 The Wilson-Fisher Fixed Point 

A more interesting case is when the solution is non-trivial. We only find one 
example of a non-trivial point, which we deem to be an approximation to the 
Wilson-Fisher fixed point. To find this fixed point we need to rely on numerical 
methods as outlined in appendix A. We will impose ( p.23| ) as a condition at the 



origin and force non-trivial behaviour by imposing ( |2.21| ) a.s (f)"^ oo . The results 
of these calculations are shown in figure (|2.1|) for various values of A^. 
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2.3 Critical Exponents at Leading Order 



The next step in the program is to calculate the critical exponents for v and lo 
corresponding to the fixed points. To calculate the critical exponents it is neces- 
sary to linearize about this fixed point potential, We will write V"(0^, t) = 
+ 51^(0^), where 51^(0^, t) is given by £:e'*'*i;(0^), with £ <^ 1, and expand 
to first order in e. Upon doing this we have 

(A - 3)t;(0 ) + 1; (0 )0 = tt^-^ttTTTTtxT^^ ) + 



(2 + 2^' (02))3/2 ^ ' (2 + 2\/'(02) + 402^// ^02)) 3/2 

(2.25) 



2.3.1 Critical Exponents for the Gaussian Fixed Point 

For the Gaussian fixed point we will expect to find perturbations with eigenvalues 
given by the classical scaling dimension of the coupling (cf Introduction), eg an 
operator corresponding to the mass with dimension 2, an operator corresponding 
to the 0^ interaction with dimension 1 etc. In the case of the Gaussian fixed 
point equation ( |2.25| ) becomes. 



(A - 3)t;(02) + t;'(02)02 = + ^ ) (2.26) 

The general solution to this equation is 

A[,_3]F(v^02)[^/,] + fi2-V2+^/^(0^)f,-_\_,^,_^/,jF(--i^)o (2.27) 
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where [ai,a2,...,ai]F{x)[b-i^^b2,...,bj] represents Barnes' extended hypergeometric func- 
tion, defined by, 

V^^"^=l r(a^) ) (C) r,o\ 

[a„...,a^]F{x)ib,,...,b,] = 2. -TlTj — r{b,,+k) \ ^ (2.28) 

If any of the = —k, for k a non-positive integer, then the series stops after 
k terms. If we wish the solutions to be bounded by polynomials then we are 
forced to set the coefficient B in the above to zero and chose A — 3 = fc, for some 
negative integer k. This then forces A = 2,l,0,— 1,..., and we have the expected 
results. We see that the at the Gaussian fixed point the scaling dimensions are 
given by the canonical dimensions [EO], that r] = and u = 1/2. 



2.3.2 Critical Exponents at the Wilson-Fisher Fixed Point 

To calculate the exponents at the non-trivial fixed point we will need to think 
about the boundary conditions. By linearity of the perturbation we can choose 
f (0) = 1. We will again ensure that the perturbation exists at the origin by 
imposing the equation as a boundary condition at 0^ = 0. That is, 

(A - 3)f (0) = rTT^^'(O) (2.29) 

^ ' ^ ' (2 + 2l^'(0))3/2 ^ ' ^ ' 

The case = is treated slightly differently. If we were to impose ( p.29|) then we 
either have A = 3 or f (0) = 0. We discard the case of A = 3 as this corresponds 
to the uninteresting case of the vacuum energy operator. We must therefore have 
f (0) = 0. The other boundary condition again comes from linearity, which we 
take to be f'(0) = 1. 

We can also consider the solutions for large 0^. This time we will see that the 
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perturbation will be a linear combination of [cf^Y"^ and exp(|[3OAf]^/^(0^)'^/^). 
Once more we will enforce a coefficient of zero on the latter perturbation, as we 
require the perturbations to grow no faster than a power in 0^. The imposition 
of this condition will ensure that the solutions satisfying (0^)^"'^ as ^ oo form 
an isolated one parameter set. 

Upon imposing these two conditions we will have a second order eigenvalue prob- 
lem for A with three boundary conditions. Therefore, we expect to see a discrete 
number of solutions, which is indeed the case. We find only one positive eigen- 
value, which yields v through v = 1/A. The least negative eigenvalue yields 
the first correction to scaling through uj = —A. The results are summarized in 
table n. 



We see that there is quite an impressive agreement between the results gained 
from other methods and those gained at this simplistic level of approximation. 
At first there is a gradual decrease in the accuracy of the approximation for u 
and a slight improvement in uj, as N increases. This has been noticed before by 
several authors and should therefore not be entirely surprising The results 
are better in the large regime, but this should be expected as the local potential 
approximation effectively becomes exact at = oo. We will discuss the special 
case of A^ = oo in a later chapter. The results compare well with those obtained 



by others using different forms of cutoff ^ Q ^ and those obtained using a 



sharp cutoff pO| , |38| , |36[] . Using a different form of smooth cutoff will not generally 
preserve a re-parameterization invariance of the equations and so other authors 
are forced to appeal to some 'heuristic' argument to set a value for t]. 
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N 


V 







0.5961 


0.5880" 
0.5880(15)^ 
0.592(3)^ 


0.6175 


0.80(4)° 
0.82(4)^ 


1 


0.6604 


0.6300(15)° 
0.6310(15)^ 
0.6305(15)^ 


0.6285 


0.79(3)° 
0.81(4)^ 


2 


0.7253 


0.6695(20)° 
0.671(5)^ 
0.672(7)" 


0.6621 


0.78(25)° 
0.80(4)° 


3 


0.7811 


0.705(3)° 
0.710(7)^ 
0.715(20)^= 


0.7068 


0.78(2)° 
0.79(4)^ 


4 


0.8240 


0.7479(90)" 


0.7510 




10 


0.9380 


0.877'^ 


0.8910 


0.78^ 


20 


0.9577 


0.942'^ 


0.9469 


0.89^ 


100 


0.9938 


0.989^ 


0.9911 


0.98^ 



Table 2.1: Critical exponents of the three-dimensional theory for various values 
of iV. For comparison we list results obtained with other methods: 

a) From summed perturbation series in fixed dimension 3 at six-loop order p2 |. 

b) From the e-expansion at order e' [22 



c) From lattice calculations ^ . 

d) From the 1 /A^-expansion at order 1/A^^ ||56|| . 

e) From a recent lattice study |^ - 



f) From the 1/iV expansion at order 1/A^ ||58[| . 
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0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2.2: The perturbation corresponding to v at leading order in the derivative 
expansion, for three dimensional 0{N) symmetric field theory, for = 1,2,3,4. 
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Figure 2.4: The perturbations corresponding to v and u at leading order in 
the derivative expansion for three dimensional 0{N) symmetric field theory, for 
= 0. The normalization condition is now given by f'(0) = 1. 
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Chapter 3 

The Results at Second Order 



In the previous chapter we showed how to obtain the results coming from the 
leading order of the approximation. In this chapter we will report the results 
of taking the approximation to the next order and including the effects of the 
terms involving two derivatives on the right hand side of the flow equation. This 
will lead to a vast increase in the complexity of the equations and the boundary 
conditions supplied to them. We will again specialize to the case D = 3. 



3.1 The Second Order Equations 



Our starting point will again be equation ( p.5|) , reproduced here, 



(3.1] 



I roc 

-— tr / dqq 
In Jo 



D-i( q dC{q' 



C{q^) dq 



2-7] 



6"'' + C, 



1 -1 



iq, -qj 



We will again take r[0] to be expanded as in equation (|2.7| ). We will also force 
the coefficient functions to be non-singular for all 0^ at fixed points, and that the 
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perturbations about these fixed points grow no faster tlian a power of 0^. The 
equations at the second order in the expansion are calculated by substituting 



expression (|2.7| ) into equation and dropping all terms with more than two 
derivatives from the right hand side. We will first consider what happens to the 
coefficient functions multiplying terms with more than two derivatives in them. 
We see that a general term will satisfy, 

t) + {(p'H'{(P\ t) +n^H{^\ t))+neH{^', t) - 3H{(P^, t) = (3.2) 

where denotes the number of 0"'s in the term multiplying H divided by 2 and 
hq is the number of derivatives in the term multiplying H. We can no longer 
force that 7] be zero at fixed points, and this leads to the following expression for 
H at fixed points (setting D=3), 

(i+'))"^i,+"a-3 

H oc ~v (3.3) 

Hence, as > 2 and rig > 2, we see that H will be singular at the origin, under 
mild assumptions about t], eg rj > 0. We are therefore forced to set all the terms 
with more than two derivatives in them equal to zero. Hence at the second order 
of the derivative expansion we parameterize F as follows. 



r[0^] = J rf^x I v(0^ t) + -K{<p\ t) {d,rr + 2^(r , t) {^""d.rY] (3.4) 

A little bit of thought shows that this conclusion holds at all orders in the ap- 
proximation: if we substitute an expression into the equation that is of higher 
order than the order we are working at, then we are forced to set the higher order 
terms to zero BOl. 



The next step is to compute the inverse in equation ( p.l|) . This is not as straight- 
forward as in the leading order case and a more involved technique is required. 
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The steps that are required will depend upon the value of A^: the case of = 1 
is slightly different than the more general case. We discuss the more general case 
first before briefly reporting the results at = 1. 



We will regard [5"^ + C.d'^T /5<p''^'^^^^-^ 



as a differential operator: 



5"^ + C. 



5^T 



iq, -qj 



J d^x Q"^ 



where 



Q 



ab 



-iq.x 



5"'' + C. 



^jq.x 



(3.5) 
(3.6) 



Q""^ is a function of q and 0(x) and its derivatives evaluated at x. To calculate 



Q'^" we will regard C and as differential operators. Noting that we are 

performing a derivative expansion, we should expect Q*^^ ~ {5^^ + z/"^)~^ + ■ ■ 
where the expression obtained by dropping all terms containing differentials 
of 0" from C{q^)e~^'^^^J^-fTe'^'^■^. Hence noting that 



Qac (^^cb ^ ^ch\^ ^ l^^ac ^ ^ac^, J^^cfe ^ j^cdylQdb ^ ^^cd ^ ^cdyl^deQeb]^ ^3 y^. 



and that from (13.61) we have 



qab ^ ^ab _ g-*q.X(^(_n) 



6^r 



(3.8) 



we must have satisfying the following expression, 



Qab ^ ^^ab^^abyl^^^ac^^acyl I ^cdQcb _ g-»q.x^^_QN " - »q.XQCt6 . ^g g^ 



iq.X/^db 



The derivative expansion is performed by iterating expression ( p.9|) to the re- 
quired order, here up to two derivatives, starting with Q"-'' = 1/(6°''' + z/"^*), and 
remembering that C(— □) = — □. To do this we must calculate T^jTb, which is 
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found, after a long but straightforward calculation, to be. 



-2{d,r)<l>''K'd, - 4{d,r)<P''{<f)%<P''yK' 

+2r(p\d^(P''fK' 
-5-\d,rfZ - 2r{d,cp'')Zd^ 
-2^y\d^^^fZ' - 5"^(0"n0'=)Z 

-2(j)''(j)\(f)''D(j)^)Z' - 6''\<f)%(j)^fZ' 
-2r{d,<p'Wd,<j)'')'Z' - 20V(0^9^0'=)Z'9^ 
-2r<p\(p'd,(P'fZ' (3.10) 



This expression is then used in ( |3.9| ) to iterate Q""^ to two derivatives. This is an 
extremely long calculation and was performed using the symbolic manipulation 



package Form |6^. The only remaining step is to calculate the trace over the 
spin indices. This will require the calculation of (5"* + z/°*)~^, and the derivative 
of its inverses. This is easily done by noting that if + u"-^ = A + B(j)°'(f)'' then, 

sab Tz> 

and also by noting that if we have a matrix A{x) then = —A^^^A^^. 
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This calculation then yields the following equation for V, near fixed points, 

+ (1 + r;)0V' - 3y = - (1 - 77/4) X (3.12) 
A^- 1 1 



we refer to the above equation as the V equation. Similarly the equations coming 
from the (9^0")^ and parts of the action will be referred to as the K 

and Z equations respectively. The equations for K and Z are a lot longer and 
are relegated to an appendix. We will perform one further transformation on the 
equations. To enable an easier comparison with the N — 00 equation we will 
scale the equations as follows, 

4) ^ N4) (3.13) 
V = NV (3.14) 
K ^ k (3.15) 
Z = Z/N (3.16) 

This then yields the following equation for V, 

^F+(l + r;)0V-3l^=-(l-77/4)x (3.17) 

1 - l/N 



where we have dropped the tildes on V , K and Z. 
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3.2 The Fixed Points Solutions 



To find the fixed points we will again need to consider the boundary conditions 
supplied to the equations. We will again ensure that the solutions exist at the 
origin by imposing the equation at the origin. This will provide three conditions. 
We can also perform the asymptotic analysis as in the previous chapter. Once 
more we find that either the solution is the trivial Gaussian fixed point (77 = 
0, = N/3\/2, K = 1, Z = 0), that the solutions are only defined for 0^ < 0^, or 
that the solutions behave as follows as 0^ —>■ 00, 

V{(P^) ~ A, {(f)'')TT^ + . . . (3.18) 
^ Ak {(f)^)TT^ + ■ ■ ■ (3.19) 

/ 9x / n. -(1 + 2'?) 

Z{(P^) ~ A, (02)^+^ + ... (3.20) 



We will be most interested in the non-trivial solutions defined for all 0^. We can 
also study the perturbations about these solutions. Forcing the perturbations to 
grow no faster than a power, we see that the solutions satisfying (|3.18|) to (|3.2UD 



will form a discrete set, with the solutions parameterized by the values of Ay, Ak 
and Az- 

At the moment the value of rj as been left undetermined and it appears that 
we have a free parameter. However we can use the scaling symmetry to impose 
an extra condition and hence fix this extra parameter. We will take K{0) = 1, 
with other possible solutions being reached by using the re-parameterization in- 
variance. We now have a seven parameter set with seven conditions imposed 
and hence expect at most a discrete set of solutions. Again we only find two:- 
the Gaussian point mentioned above and an approximation to the Wilson-Fisher 
fixed point. The results for 77 are summarized in table ( p.l[ ) and the results for 
the fixed point solutions shown in figures ( p.l|) to (|3.3|) . 
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Figure 3.1: The Legendre effective potential at second order in the deriva- 
tive expansion, for three dimensional 0{N) symmetric field theory, for = 
2,3,4,10,20,100. 
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Figure 3.2: The K component of the wave function renormahzation of the Leg- 
endre effective action at second order in the derivative expansion, for three di- 
mensional 0{N) symmetric field theory, for N = 2,3, 4, 10, 20, 100. 
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Figure 3.3: The Z component of the wave function renormahzation of the Legen- 
dre effective action at second order in the derivative expansion, for three dimen- 
sional 0{N) symmetric field theory, for N = 2,3, 4, 10, 20, 100. 
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3.3 The Critical Exponents at Second Order 



To calculate the critical exponents we will again need to consider the perturba- 
tions about the fixed point solutions. That is we will write, 

= l/*(02)+£e^*t;(02) (3.21) 

= K*(0')+£e^*A;(0') (3.22) 
Z(02,t) = Z*{(^^) + 6Z{(P^,t) 

= Z*(02)+£e^*^(02) (3.23) 

where y*(0^), K*{(f)^), Z*{(f)^) are the fixed point solutions calculated above, and 
expand the equations to first order in e. 

To find the perturbations we will again need to consider the boundary conditions. 
We will want the perturbations to exist for all 0^. If we insist that v,k and z 
grow no faster that a power at large 0^, then asymptotic analysis will show that 
V, k and z will grow according to their scaling dimension, 

v{(t>'') ~ a,((/)2)Ti?-^ + ... (3.24) 
fc((/.2) ~ a,(02)T^-A^... ^325) 

„ „ l + 2rj , 

z{(i>^) - a^(02)"^"^ + --- (3.26) 

The imposition of power law growth will again force a^, and to form an 
isolated three parameter set. The other three boundary conditions come from 
forcing the equations to hold at the origin. Using linearity to set ^;(0) = 1, we 
will have a seven parameter set with seven boundary conditions imposed. Hence 
we expect to find at most a discrete set of solutions, which is what is found. As 
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before we find just one positive eigenvalue, which yields v through v = 1/A, and 
determine the least negative eigenvalue, which gives the first correction to scaling 
exponent through lo = —A. These values are shown in table (|3.1|). 



It is important to recognize that we will also find other solutions of the equations 
which do not correspond to critical indices These solutions are known as re- 
dundant perturbations and the eigenvalue corresponding to the solution depends 
on the exact form of the renormalization group chosen, but as their eigenvalue 
depends on the exact form of the equation chosen these perturbations can not 
be physical. The redundant perturbation reflect invariances of the equations. 
In general if we change variables = + £$"[0] in ( p.l|) (with J.0 re- 

placed by J.(p), then this induces a change in the effective action of ST = F.^ 
with F[(j)] = eexp (— iy[J]) ^[6/6J] exp (W^[J]) and a change in the cutoff func- 



tional fS^, |6^. A general choice of F that will leave (|2.7|) invariant at this order 
is, 

m = {fi<p\^))r + ax>^d^r{x)} (3.27) 

for any function / and any constant a. A redundant perturbation must then 
satisfy, 

V oc 2/0V'-3a\/ (3.28) 
k oc 2f(t)'^K' + 2fK ~aK (3.29) 
z oc 2/02^' + 4/'02^ + 4/Z - + 2/'ir (3.30) 



In fact we expect and find only one redundant perturbation corresponding to the 
re-parameterization invariance ( p.8| ). This has eigenvalue A = 0, / = 5/2q; and 
a ^ 0, yielding {v, k, z) = (-50 V + 3^, -5(p'^K' - AK, -h<p'^Z' - 9Z), where we 
have set a = —1. An obvious question is why no redundant perturbations are 
found at the leading order. The answer is simple: the choice of i^' = 1 breaks the 
re-parameterization invariance, and we should not expect to see any redundant 
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perturbation. We should also not expect to see any perturbation corresponding 
to 0-translation invariance as the choice of V, K and Z being functions of cfP' 
breaks this invariance. 



The redundant perturbation also proves to be a convenient check on the numerical 
accuracy of the equations for the perturbations. Once we know the fixed point 
solutions we automatically know a solution to the perturbation equation, which 
should have eigenvalue A = 0. We can then try to find our redundant solution 
using numerical methods. The degree to which the known solution (from the 
fixed points), and the solution found numerically agree, gives a bound on the 
accuracy of the critical exponents. Using numerical methods, we find that the 
value for the eigenvalue corresponding to the redundant perturbation typically 
lies in the range 0.0003 (for = 2) to 0.00005 (for N = 100), indicating that 
we should only trust any numeric results to three figures. The graphs of the two 
types of solutions also agree to a high degree of accuracy. 

3.4 The case of iV = 1 

As mentioned above the case of = 1 is slightly different to the more general 
case oi N ^ 1. This should not be entirely surprising as there is now no internal 
symmetry at all (other than the discrete Z2 symmetry). The derivative expansion 
now becomes, 

r[0] = j d^x [V{ct?) + \ {K{ct?) + 0^Z(02)) (9,0)^} (3.31) 

We see that we should no longer consider the Z and K components of wave 
function renormalization separately, but should consider instead a single function 
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Figure 3.4: The results of applying the derivative expansion at second order for 
the case of iV = 1. 

K = K + (fP'Z. The derivative expansion at second order, at = 1, thus becomes, 



(3.32) 



We will no longer have three coupled second order equations but will have only 
two coupled second order equations for V and n. The results for this case can 
be found in It is an important check on our equations that if we write 

n = K + (fP'Z and set = 1, then we get the same equations as found in |]30 



The results of this expansion are included in figure (3^) for completeness. These 
come from an independent calculation. 

It is interesting to ask what does happen to the K and Z components as — 1. 



This can be done by taking the K equation and substituting K 



K 



into it. This then yields an equation which can be solved numerically using the 



known result for k and V (from [^| or our independent calculation). It is then 
found, that at = 1, Z diverges as 1/0^ at the origin. Again this should not be 
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entirely surprising: we can write Z = — K)/(fP'^ and should therefore expect 
to see some divergence at the origin, provided -ft'(O) 7^ /t(0). Numerically this 
is an important point. As ^ 1 we will expect that Z will become steeper 
and steeper at the origin before it eventually becomes divergent at = 1. This 
divergent behavior may lead to numerical instability as we approach the regime 
A^ ^ 1. In fact we start to feel the effects of this divergence as early as A^ = 2, 
where it is already found to be much harder to produce an accurate numerical 
solution than at A^ = 3, say. 



A^ 




V 


















0(9°) 




1 


0.054 


0.035 


0.62 


0.66 


0.63 


0.898 


0.63 


0.80 


2 


0.044 


0.033 - 0.04 


0.65 


0.73 


0.67 


0.38 


0.66 


0.79 


3 


0.035 


0.033 - 0.04 


0.745 


0.78 


0.71 


0.33 


0.71 


0.78 


4 


0.022 


0.025 


0.816 


0.8240 


0.75 


0.42 


0.75 




10 


0.0054 


0.025 


0.95 


0.93 


0.88 


0.82 


0.89 


0.78 


20 


0.0021 


0.013 


0.98 


0.96 


0.94 


0.93 


0.95 


0.89 


100 


0.00034 


0.003 


0.998 


0.994 


0.989 


0.988 


0.991 


0.98 



Table 3.1: Critical exponents of the three-dimensional theory for various values 
of A^. We list the results gained at second order together with those from the 
leading order and an indication of what should be expected from other methods 



as summarized in p2|, ^ ^ . 77 is predicted to be zero for all A^ at the leading 
order in the approximation. 



3.5 Discussion 

The results at second order are puzzling: for A^ = 2, 3, 4 we have reasonable 
accuracy for ri and z/, with an improvement upon the leading order results, but 
terrible predictions for tu, which is predicted to be worse than the leading order 
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results. On the other hand the predictions at = 10, 20, 100 show the converse: 
we have rj seriously underestimated, by a factor of about ten, and the results 
for V compare less favourably with other methods than those predicted by the 
leading order of the approximation. However we also see an improvement in the 
predictions for uj compared to other methods. 

These results are very puzzling as the only other comparable calculations at this 



order in the derivative expansion I^O], |6J], performed at = 1, showed an im- 



provement for all the calculated exponents when the change was made from the 
leading order to the second order of the derivative expansion. Other authors who 
have tried to go beyond the leading order of the approximation certainly do not 
report results for oj like those we have found at = 2, 3, 4 |61|]. It would be easy 
to explain these results away by calling upon numerical inaccuracy, especially at 
= 2,3,4 where the known steep behaviour of Z at the origin may have an 
effect. However, although it is true that the equations are extremely difficult to 
solve (cf appendix A) we feel that this is unlikely for several reasons. Firstly the 
results are insensitive to changes in where we impose the asymptotic boundary 
condition, the number of points in the solution mesh and other numerical factors, 
whereas we might expect there to be random errors introduced if the results were 
down to numerical roundoff or discretization error. For example, it can be seen 
that the eigenvalues seem to flow continuously with the value of N, whereas we 
would perhaps expect some random fluctuations if we were encountering severe 
numerical problems. Secondly, we already know that the known redundant per- 
turbation can be found to a reasonable degree of accuracy, indicating that we 
should have a high degree of faith in our numerical methods and in the numer- 
ical stability of the equations themselves. For these two reasons it seems rather 
unlikely that these effects are due to numerical problems. 

Quite why the results do not compare so well with those presented by others is 
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• Derivalive expansion at second order 
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• Derivative expansion at second order 
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(a) (b) 

Figure 3.5: The flow of rj and u witli tlie value of N. Figure b also shows what 
the large N expansion predicts at order 1/N. 

puzzling. The serious under-prediction of rj for large could be attributed to 
the approximation scheme being 'too biased' towards the N = oo results, which 
would predict rj = 0, at large A^. It is also interesting to note that the flow of u 

with seems to follow the large results at order 1/A^ for A^ = 3,4,5, 

The results of these two flows are shown in flgure (|3.5| ). 

Overall we must conclude that the scheme for intermediate A^ does not work as 
well as it does for A^ = 1 or A^ = oo. The N = oo result is in fact exact and is well 
understood (cf the next chapter and The reason why the encouraging results 
at A^ = 1 [^, Q are not continued for A^ 7^ 1 are not yet understood. However, it 
should be emphasized that we flnd no evidence for spurious results, such as those 
found by truncations: we flnd only the flxed points and perturbations expected. 
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(c) z{ct>^) 

Figure 3.7: The f , k and z components of the perturbation corresponding to uj 
at second order in the derivative expansion 
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Chapter 4 



Exactly Solvable Systems 



It has been well known that certain theoretical models can be exactly solved. 
That is, have all their physical parameters determined without any form of ap- 
proximation. In this section we will discuss how the derivative expansion repro- 
duces the exact results found in two well known cases:- N = (yo and N = —2. It 
is important that the derivative expansion reproduces these two exactly known 
results, for if it didn't then the approximation would be too great a mutilation 
of the theory. 



4.1 Exact Results at = oo 



It was noted a long time ago that the large limit is equivalent to the 
exactly solvable spherical model pBI. It was shown that in the large iV-limit that 
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the local potential expansion effectively becomes exactQ ie we have, 

T[^] = J d^'x^idf.cP^f + V{(P^) atAr = cx) (4.1) 

In this section we will show that the derivative expansion re-produces this result 
and give an expression for V^(0^) and for the perturbations about the fixed point. 

In the previous chapter we scaled the equations to provide a comparison with the 
N = oo results. We now justify where these scalings come from. 

We will consider a simple scalar $^ theory in D dimensions, 

Z = y"D$exp|-5[<l>] + y"ci^a;J.$| (4.2) 

We can perform the following transformation 

exp d-. l<tr ~ / exp {/ d".^ - -^<S>'^^ (4.4) 

and hence re-write the path integral as. 



Z = I^Aexp || rf^x (^iVA^) + ^trlog ^-A + + | 



(4.5) 



We see from this that we now have an explicit dependence on in the (non-local) 
action. If we wish to perform a saddle point expansion in 1/A^ then we must have 
a factor of outside the action, but none in the trlog. From this we see that 
we must hold Nu fixed, so u is of order 1/A^, and we must have A of order one 
also. Hence, counting the powers of in equation ([4.4]) and demanding that the 



^Wegner pointed out that although the other terms exist in the expansion, they do not affect 
the value of the critical exponents M 
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powers of match up on both sides, we see that we must have $^ of order A^. 



From this we see that when we come to consider the flow equations we should 
scale 0^ as 0^ = Ncf)"^. Hence, demanding that F scales as F = A^F (as follows 
from (|4.5|) , we have the following scalings, 

V ^ NV (4.6) 
K ^ k (4.7) 
Z ~ Z/N (4.8) 

where V ,K and Z are of order one in the expansion. Upon doing this the 
flow equation for V becomes (upon dropping the tildes), 

+ (1 + r/)0V' -W = -{l- r//4) X (4.9) 

1 ^ 1-1/A^ 



N^K + {2V' + 402 V" + 2^K + 02 ^'"^ (2V' + 2/^^ ^'"^ 



plus some longer equations for K and Z. If we now take the N = oo limit we 
end up with the following sets of equations, 

|y + (l + r^)0V-3y = - ^/^""^^^^.i/. (4.10) 

l-K + {l + ,)<P^K' + ,K = - ^^/'^^'^^,3/2 (4-11) 

( 2V' + 2v^ ) 



^Z + (l + r/)0V + (l + 2r/)Z = ••• (4.12) 



We immediately notice something different:- the equations for K and V have de- 
coupled from the equation for Z and the equations are now first order equations. 
This is the first simplification. 
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Now lets consider the fixed point solutions. The next simplification comes by 
looking at the local potential expansion and noting that this will require that 



K = 1. Hence we try K = 1 as a fixed point solution. Equation (|4.11| ) then 
yields, r] = 0. 

Upon setting r/ = we find that the fixed point potentials will then satisfy, 



/,V -3V = Tj. (4.13) 

{2V' + 2)'/' 



Although this is a first order equation it is still not in a particularly useful 
form for solution. If we differentiate the equation with respect to 0^ and set 
H/(02) = ■l/'(02) then we have. 

This equation can the be re-arranged as partial differential equation for 0^ in W 
as follows. 

Using as an integrating factor we then have, 

1 



dW 2W^I\2 + 2W^)3/2 
Performing the final integration then yields a final expression for 0^, 

W 



(4.16) 



02 = ^^+^!^^ \ (4.17) 

where A is a constant of integration. If we have A 7^ then it can be shown that 
0^ will reach a maximum value for some non zero value of W . ie the solution is 
not defined for all values of 0^. We are therefore required to set A = and the 
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solution then becomes, 



2v^ 



VT+wii 



w 



l + W 



This can be inverted to give the solution for W, 



W 



-1 + 



1 

V2 



1 2 



k2\2 



+ 1 



Integrating with respect to (j)'^ then gives the solution for the potential, 



V 



(4.18) 



(4.19) 



(4.20) 



It is important to stress the difference between the scaled and the unsealed results 
in the large N regime. The unsealed V tends to infinity, K tends to 1 and Z 
tends to zero as oo. However when we consider the scaled solutions we see 
a different pattern: V now tends to our solution ( [4.20| ), K still tends to one, but 
so far Z has been left undetermined. 



77 



Performing the scaling upon the Z equation and setting X = 1, = we see 
that Z must satisfy the following equation, 



fZ' + Z = 



+ 



■' I ^2 T.// 1024 4 

+ (4 Z6,2,3 - 8 ^4,2,3 + «4,2,2 + 4 ^2,2,3 " 5 ^2,2,2 + «2,2,l) 

64 2 / 2 
+ y0 (^16^5,3,2^0 - 8«5,2,3 - 11 «3,2,2 + 8«3,2,3 - 8Z3,3,2 

+8 ^3,3,l + 8^5,3,2) V^'" 

/256 27. 256 2 7 . ^ 256 2 7 • 

V~3~ ^'^'^ 3~ ^"^'^'^ ~3~ ^"^'^'^ 

416 „ ^ 256,^2x2^2 128 256 

3~ ^ ^ ^4,2,2 + ) ^ ^.3.2 + ^6,2,3 ^ «4,2,3 

128 64 64 128 

- 64 Z4^2,2 + «2,2,3 - Y «2,2,2 - y «2,2,1 + «6,0,5 

64 . 128 . 128 . 256 . 80 . 

+ «6,3,2 + 4,3,1 ^ «6,1,4 + 4,1,4 + 4,1,3 

128 . 16 . „^ . 64 . 64 . 

H ^ ^,0,5 + -g- «2,0,3 — 32 Z2,o,4 + ^2,3,0 + ^2,3,2 

128 , 112 . 16 , 128 . 32 , 

^ ^2,1,4 H — ^,1,3 + ^2,1,2 — ^,3,1 ^ 4,0,4 

256 . 128 . 256 ^ 256 2 ^ • \ ... 

- «4,0,5 - «4,3,2-^ ^ «4,3,1 ~ ^ ^.2,3 1 

^ ^8 ^2 7/ . ^8 40 2 7, . 64 . 

+ ( - / Z3,2,l + - Z Z5,2,2 - y / ^3,2,2 - y Z Z5,2,3 

40 ^ . ^8 ^ 32 64 

+ g ^ «3,2,1 + y ^ ^3,2,3 ~ y ^ ^5,2,3 



rll2 



40 ^ . ^8 

-y^^3,2,2+g. 



+ 



y 0^ ii7,2,3 + 8 Z ^3,o,3 + y ^ ^7,2,3 + ^ ^' ^3,2,3 

+ \<t?Z' 4,2,2) 



where ij,fc,m represents the following integral, 



/A 



(1 + (21/' + 402\/")?2 + (1 + 02^)^4)^ (1 + 2V'q^ + 



Of course this equation can no longer be solved analytically, as above, but the 
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Figure 4.1: The potential, V , and the Z component of the wave function renor- 
mahzation of the Legendre effective action at second order in the derivative ex- 
pansion, at = cxD. 



numerical solution is straightforward. Imposing the constraint that Z exists for 
all (jP' we get the solution shown in figure ( ^4.1| b). 

We have now calculated what the fixed point potentials look like, but what we 
are really interested in is the critical exponents. To calculate these we will again 
linearize about the fixed point solutions by writing = +£:e'^*f (0^) 

and K{(j)^,t) = K*(0^) + £e^*A;(0^), where V* and K* denote the fixed point 
solutions found above. Using these expressions in equations ([4.10|) and (^4.11|) and 
expanding to first order in e yields the following pair of differential equations. 



Xv + (j) V — 3v 
Xk + (b^k' 



1 k + 2v' 
k' 

{2V' + 2)3/2 



(4.21) 
(4.22) 
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We will consider equation ( |4.22| ) first. Re-arranging the equation we have, 

X{2V' + 2)3/2 



()2(2l^' + 2)3/2 _ I 



(4.23) 



Using expression ( [4.20|) we then have 



2. (A2 



k' _ A {2(02)2 ^ 1 ^ 

^ 02 |2(02)2 + 1 + 202^^(02)2 + 11=^/' 

This can easily integrated to give the solution for fc(0'), 



3/2 



(4.24) 



(4.25) 



where Ai is a constant of integration. 

We now consider the v equation. This can be re-written as follows, 

(A-3)(2 + 2n^/2. 1 k 

02(2 + 2V')^/' - 1 2 02(2 + 2\/')^/' - 1 

Using an integrating factor of (exp(2 sinh"^ 02) — 2)'^/2^3/2^ -|-]-^jg equation becomes. 



d 



t;(02) (exp(2 sinh^i 0^) - 2)^/2-3/2 



lA; (exp(2sinh-^ 02) - 2)^2-3/2 
2 



\2 + 2^')^/' - 1 



(4.27) 



Using the known results for V{(fP') and fc(0') this equation can be integrated to 
give. 



+A2 



- 2 



-A/2 



g2sinh-^^_2)=^/2-^/' (4.28) 



From equation ( [4.25| ) we note that for A;(0) to be real we need A = —2m, where 



m 
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is an integer, or have Ai = 0. However, we also note that to avoid a singularity 
near the origin we must have m positive or zero ( as exp(2 sinh~^ 0^) — 2 < for 
0^ < sinh(log v^))- On the other hand, considering equation ( [4.28|) we see that 



we must also have 3/2 — A/2 to be a positive integer, or have A2 = 0. We see 
that we must have two mutually exclusive sets of solutions:- 

• a set of solutions with eigenvalues A = 3 — 2m, for a positive integer m, and 
Ai = 0. These solutions have A = 3, 1, — 1, . . . and satisfy 

^(02) = ^2(e2sinh-^^_2)'^'"'^' (4.29) 
A;(02) = (4.30) 



• a set of solutions with A2 = and A = —2m, for a positive integer m = 
0, 1, 2, 3. These solutions have A = 0, —2, —4, . . . and satisfy, 

^ ^^_L (^2e3-h-i02 _ ^ 3e-^'°'^"'*') 

^(^g2sinh-l<^2_2)-^/2 (4.31) 

fc(02) = Ai(e2^''^^-^'^'-2)"^^' (4.32) 

The first set of results reproduce those first found by Wegner and Houghton |]4l[| , 
whilst the second set correspond to corrections to scaling coming from the (9^0")^ 
part of the action. We plot some of these functions for various values of A in 
figure (|4.2|) . It should be noted that all the solutions meet at the point (0, 1/2^/2). 
This behaviour is related to the singularity structure of the equations and is not 
of any particular importance. 
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(a) v{(lP) for the first type of solutfon (b) v{4>^) for the second type of solution 




(c) for the second type of solution 

Figure 4.2: The v and k components of the relevant, leading irrelevant and sub- 
leading irrelevant operators at = oo. 
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4.2 Exact Results at iV = -2 



The fact that the case of = —2 is also an exactly solvable case was noted long 
ago |6^. The initial hope was to be able to combine the known results at 
N = oo with the new results at = —2, and use Pade approximants to gain 
information about the physically interesting case of = 0,1,2,3 O, IBSI, IBD 



It was found that when = —2 the exponents reduce to those of the Gaussian 
model, ?7 = 0, z/ = 1/2. However there is one large difference: the A^ = —2 
case describes the critical Gaussian model, as opposed to the tricritical Gaussian 



model fTy]. Therefore, we only expect to find one relevant eigenvalue instead of 



the two relevant eigenvalues that we found in section 2.3.1 



This phenomena can be understood as follows. As pointed out by Fisher ||68|| , 
at A^ = —2 only the quadratic part of the action plays a role, ie the physics is 
totally independent of any quartic coupling etc. At the tricritical fixed point we 
see two relevant operators corresponding to a mass term with dimension two and 
a (j)^ term with dimension one. However, as we are now reduced to only having 
the quadratic part of the action playing a role, we only can expect to find one 
relevant operator. 

We will show how the derivative expansion reproduces these results at the lead- 
ing order of the approximation before briefly discussing what happens when we 
consider the second order terms. 

Our starting point will be the equation for the potential at leading order, 



A^- 1 



A V -W = -{ y ^ + ^ — — \ (4.33) 



It turns out that we need only consider what happens at the origin. If we differ- 
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entiate the equation and set 0^ = we get 

N + 2 



(2 + 2y'(0))3/2 



y"(0) - 2^' (0) (4.34) 



We see that at = —2 this equation simphfies to — 2V(0) = 0, so we must 
have V(0) = 0. Considering equation ( |4.33| ) at 0^ = we see that we must 



then have ^(0) = — -\/2/3. From this we see an interesting departure of the 
behaviour at the origin from the cases = 1,2,3,..., etc. Beforehand the 
boundary condition at the origin used to leave V"(0) and V'{Q) as free parameters, 
constrained by the given boundary condition. However we now have 1^(0) and 
V(0) exactly determined, with V^"(0) decoupling and becoming a free parameter. 
The question is whether we can find a non-trivial solution: after all if V(0) = 
then series methods show that we only reproduce the tricritical Gaussian solution, 
V^cjP') = N/{?)\/2). All analytic attempts to find such a solution failed, but if 
we use ^(0) = — v^/3 and V(0) = as boundary conditions then numerical 
methods quickly find a non-trivial solution with V""(0) 7^ satisfying the required 
asymptotic conditions at large 0^ (cf equation (|2.21| )). The results of this search 
are shown in figure (^]3|a). 

We now consider what happens to f (0^). We consider the linearized equation, 

i>^-^> + 21/03/2 + (2 + 2V' + A<P^V"fl- ^^-''^ 

If we differentiate this equation and set 0^ = we have, 

(iV + 2K 3'^ + ^'";!::+. '(A -2) (4,36) 



(2^/ + 2)3/2 (2\/' + 2)5/2 

Hence we see that if we set A^ = —2 then the terms involving V" and v" decouple 
and we are left with, 

v{X - 2) = (4.37) 
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(a) The fixed point potential, ¥{4)^) (b) The perturbation v{(j)^) for A = 2. 

Figure 4.3: The fixed point solution for tlie critical-Gaussian fixed point, together 
with the only found perturbation, at iV = —2. 

Therefore, we either have A = 2 or we must take f'(0) = 0. Hence we see that we 
find the first eigenvalue without re-course to any analytic or numerical method, 
and so v = 1/2. Again we check numerically that we find a solution that isn't 
just the tricritical Gaussian solution. This is shown in figure ( [4.3| 3). If we take 
v'{0) = 0, then setting 0^ = in the linearized equation at = —2, shows that 
f (0) = 0. We now have a second order equation with two boundary conditions 
set at the origin and hence we expect to find a unique result. By inspection we 
see that this solution must be vid)"^) = 0. 



As pointed out by Fisher [£8[ we should also expect to find exact solutions for 
A^ = —4, —6 — 8, . . .. We find this to be the case, eg consider A^ = —4, following 
a similar method to above, except this time differentiating twice, yields, 

^^ (.-i)(,^-,y(0) ^^^^^ 
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Prom this we see that we must have A = 2 or A = 1 or have f'(0) = 0. As before 
the case of f'(0) = yields no further results. This result should be expected 
because, as Fisher pointed out ||6^, at = —4 interactions up to degree 0^ now 
play a role. The program can be carried out for any negative even integer to yield 
similar results, so for = —6 we have A = 2, 1, and interactions up to degree 
0^ become involved. 

When we come to perform the analysis at the second order the equations are 
vastly more complicated. Due to the fact that even modern symbolic manipula- 
tion packages cannot handle the equations efficiently we report no further results. 
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Appendix A 
Numerical methods 



In this section we will outline the methods used to solve our equations. We begin 
with the numerical methods used before discussing particular problems that we 
encountered. 



A.l Numerical Methods 

The equations form a two point boundary value problem. There are two such 
methods of solving such equations: shooting |7^, and relaxation [|77|. Both 
methods were used at some point in the research. 

A. 1.1 Shooting 

Suppose we have a set N of differential equations, 

dyi{x) 



dx 



gi{x,yi,...,yN) i = l,2,...,N (A.l) 
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with two boundary conditions. At the first boundary, x = Xi the solution is 
supposed to satisfy, 



Bj{xi,yi,...,yN) = j = l,...,ni (A.2) 
whilst at the second, x = X2 it is supposed to satisfy, 

Bj{x2,yu...,yN) = j = l,...,n2 (A.3) 

with rii + n2 = N. 

In the shooting method we choose a set of values for all the dependent values at 
the first boundary. These values are chosen to be consistent with the boundary 
condition (|A.2|) , but are otherwise free to depend upon a set of free parameters, 
the values of which we guess. Using the values of the dependent variables as ini- 
tial conditions we integrate equation ( [A.l| ) out to the second boundary, using a 



standard ordinary differential equation integrator. Of course we will usually find 
that there is a discrepancy in how well the values at the second boundary fit the 
boundary condition there. How well the values at the second boundary condition 
fit the boundary condition depend upon the initial choice of the dependent vari- 
ables, so we now have a problem of finding what initial values of the dependent 
variables fit the second boundary condition the best. Hence, the problem has 
been reduced to a multi-dimensional root finding problem which can be solved 
using any of the standard methods of numerical analysis (eg Newton-Raphson 
etc). 
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A. 1.2 Relaxation 



In the relaxation method the equations are replaced by a finite difference mesh 
that spans the range that wc arc interested in. An initial guess is then supplied: 
this guess need not satisfy the equation or even the the boundary conditions. An 
iteration process known as relaxation is then performed to bring the equations 
into closer agreement with the finite difference equations and the boundary con- 
ditions. This again boils down to solving a multi-dimensional set of non-linear 
equations using some numerical method. 

Relaxation is particularly suited to equations and boundary conditions involving 
complicated expressions which can't be solved in closed form, ie there is no need 
to explicitly write the equation in the form dy/dx — f{x,y). It is also the best 
method when we need to find solutions that depend upon some parameter, such 
as the value of N in our 0{N) symmetry. Once we have found a solution for 
one value of we can use this solution as an initial guess for a close value of A^. 
Given the nature of our problem it is not surprising that relaxation turns out to 
be the principal method used during the research. 



A. 2 The Boundary Conditions 

Before discussing the solutions of the equations it is necessary to discuss the 
equations and the boundary conditions supplied to them. 

We wish to formulate the equations as a set of non-linear, coupled second order 
differential equations. Noting that the equation relating to the Z part of the 
action contains powers of V" (cf appendix B) it is necessary to differentiate the 

V equation and find an expression for V" in terms of V.V ,V" , K, K' , Z and 
Z'. Similarly, when we come to look at the boundary condition relating to the 
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Z equation we see that it contains powers of V" . Again we differentiated the 
V equation to find an expression for V^"(0) in terms of V", V\ K and Z. Similar 
transformations where made when we determined the equations for the pertur- 
bations. 

The other set of boundary conditions seem to be more problematic. We need the 
solutions to exist for all cfP' . We know that if the solutions exist for all 0^ then 
we can work out how the solutions behave at large 0^. This is what provides the 
second boundary condition. We choose a value of cfP' such that the asymptotic 
expressions (shown in appendix B) become a good enough approximation to the 
solutions. We then use the asymptotic expressions for V, K and Z to provide 
boundary conditions for V, K, K', Z and Z'. The value of 0^ chosen, which 
we call 0^sy, must be chosen large enough to ensure that the solutions satisfy the 
asymptotic conditions sufficiently well, but not so large that numerical instability 
prevents us from obtaining a solution. There is a certain amount of trial and error 
in deciding where to set 0^^^: we cannot really decide where to set it until we 
know something about the solution. The solutions were also checked to be stable 
against reasonable changes in the value of 0^y. 

In fact it is necessary to develop the asymptotic expressions to beyond the leading 
order in order to find a solution, eg consider equation ( |2.2(J| ), if we just use 
the leading order results, where the solutions behave according to their scaling 
dimension, then we see that, as we set (pl^^ at a finite value, then V"{(j)'^) would 
be left undetermined. Hence, we are forced to expand the asymptotic expressions 
to beyond leading order. A similar result holds for K{(f)'^). However, when we 
consider Z{(j)'^), we see that we need to take the asymptotic expressions to next 
to next to leading order. This is because when we calculate the correction to 
the leading order asymptotic behaviour it does not involve Z" , hence leaving Z" 
undetermined. In order to ensure that Z" gets determined, we need to extend 
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the asymptotic expression to one further order. 

When we come to consider the perturbed equations it is sufficient to use the 
leading order asymptotic expression. As a check we calculated the corrections to 
the leading order behaviour and also used these. This made no difference to the 
values quoted. 



A. 3 The Solution of the Equations 

Having determined the boundary conditions we can discuss how the numeri- 
cal solution of the equations was performed. The solution of the equations is 



complicated by the fact that they are stiff |7§, [T^]. This arises because small 
perturbations from the true solutions lead to a singular behaviour. Relaxation is 
particularly adept at handling stiff problems. Before discussing these issues we 
briefly discuss the problems posed by the integrals in the equations. 



A. 3.1 The Integrals Ia^b,c 

The integrals Ia,b,c play an important part in the solution. Defined by, 

^ /-oo g2a 

"'^'"~Jo (1 + {2V' + A<P^V")q' + {K + <p^Z)q^f (1 + 2V'q^ + Kq^f 

it is possible to find a general analytic result for them. This proves not to be 
a sensible thing to do. The integrals are generally small numbers, but in the 
analytic expressions these small numbers arise from the cancellation of two large 
numbers. Roundoff errors start to play a large part in determining their value 
using the analytic expression. Even performing the calculations using higher 
precision Fortran routines could not provide accuracy beyond four significant 
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figures. This level of accuracy was totally insufficient and so numerical methods 
were used. 

We used an adaptive integrator from the NAG libraries. The problem was split 
into two parts: a part from zero up to some finite value, and the remaining bit. To 
avoid roundoff errors we used an expression arising from asymptotic analysis to 
calculate the second part of the integral, and only used the numerical routines in 
the part of the real line from the origin. The point were the crossover took place 
was set to be as small as possible (to avoid numeric round off in the numerical 
integration), whilst still allowing the asymptotic expression to be accurate. 



A. 3. 2 The Z equation at = oo 

The equation that determines Z at = cx), equation ( [4.21|) is a first order non- 
linear equation. Being first order we only have one boundary condition, which 
is that the solution exists for all 0^ and has the required asymptotic behaviour. 
For this problem we used the shooting method. An initial value of Z was used 
to integrate the equation out towards the second boundary using an eight point 
Runge-Kutta integrator. In general the initial value of Z would be wrong: if it 
was too large then the solution would tend to infinity at some finite 0^; if it was 
too small the solution would tend to minus infinity at some finite value of 0^. 
Using a binary chop we could home into the correct value of 2'(0) and find the 
solution that was defined for all 0^. 

It was particularly difficult to get past the point 0^ = l/2-\/2. This is because 
at this point the right hand side of the equations has a 0/0 type behaviour, ie 
if Z' = n(02, Z)/c/(02, Z), then as 0^ 1/2/v^ we see that n(02, Z) ^ and 
d{(f)'^,Z) 0, for the true solution. There are general methods of dealing with 



such singular points |7^, but our equations were sufficiently well behaved to avoid 
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needing to use them. 



A. 3. 3 The Numerical Solution of the Flow Equations for 
General A^^ 

The numerical solution of the leading order equations proved to be easy. It was 
only necessary to provide an initial guess and the relaxation routines quickly 
found the solution. 

When we come to consider the second order equations we see that there is vast 
increase in the complexity of the equations. The numerical solution of the second 
order equation is made difficult by the large number of integrals that need to be 
calculated during each iteration. This places a huge computational overhead on 
the numerical calculations. It was soon noted that the computational power re- 
quired could not be provided by normal machines. Parallel Fortran code, making 
use of MPI [|T8[ , was developed to run on a sixteen node IBM SP2 system. The 
use of a parallel machine considerably speeded the calculation up, with run times 
typically between ten minutes to one hour (before developing parallel code the 
programs could take up to one week to run). Using the N = oo solutions as an 
initial guess allowed us to solve the general N cases using relaxation. 

The shooting method was tried, but this proved to be unsuitable. This is because 
to use shooting effectively we need to be able to reach the other boundary or have 
a particularly simple set of equations (eg the N = oo equations) . The singularity 
structure of the equations makes it impossible to shoot from one side to the other 
and hence makes the more complicated problem unsuited to shooting. 

The main difficultly in solving the second order equation was removing the insta- 
bilities in Z near the origin. These manifested themselves in the form of sharp 
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spikes in Z' near 0^ = 0. These were removed by writing ^(0^) = 

where s((/>^) is a known scahng function that makes as flat as possible. This 

greatly speeded up the calculations and removed the spikes in the solutions. 

When we came to consider the equations for the perturbations it was realized 
that we could calculate the parts that relied upon the fixed point solutions and 
store them in a file, thus saving a large computational overhead. It was important 
that these functions, which are the functions that multiplied v', v", k etc in each 
of the linearized equations, are well determined. The size of the z equation was 
such that small errors in the fixed point solutions made it difficult to determine 
these functions accurately, for large 0^ . This problem was solved by noting that 
we could determine these functions by using the known asymptotic expressions 
for y, K and Z and then matching them onto the ones calculated from the fixed 
point solutions. It was checked that the eigenvalues were stable against reasonable 
perturbations in where 0^ was set and the number of points in the mesh. 
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Appendix B 

The Equations at Second Order 



In this appendix we show some of the expressions thought to too long to appear 
in the main body of the text. We include the equations for K and Z as well as 
the expressions for the asymptotics to next to leading order. 



B.l The Second Order Equations 

In the the expression below Ia,b,c represents the following integral, 



roc 

^a,b,c ,^ ^ /mri , a i9t ^.9 , / , 19 r7\ ^A\b 



These integrals are not evaluated as the equations become very large. At present 
there are about three hundred terms in the equations. Evaluating the integrals 
raises this to between twenty five and thirty thousand terms. 
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The K equation is given by 



{1 +riWK' + r]K 



(4-r;) 
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The longer Z equation is given by, 



(1 + 77)0' Z' + (l + 2ri)Z ^ 
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(32V"'K 8K" 32 o 8 



4,0,3 I 3 ^ N 3 ; 



+/450 - 384— 384 — ^ - 512 — ^ - 512 — 

512 V'"'^ K {(p'^f 512 16 K'^ 32 K' Z 

3 N 3 TV ^ y iv~ ^ y ~/\r 

32 K' Z'ct)'^ 16 Z2 32 Z Z' 16 Z'"^ U'^f 

H ^ H ^ ^ H ^^-^ 

3 A^ 3Ar 3 A^ 3 A^ 



/ K'^ 56K'Z bQK'Z'(f 28 56 Z Z>2 

+ /4, 3 ,0 (^4 — + - ^ + - + - - + - 

28 (0^)^ 16 2 
3 A^ 3 



/ 80 F;;^^ 28if'^ 80 2^ 28 2^ 
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, 8 if'^ K'Z 5 ^2 4 



3,0,4 I 3 ^ 3 



/ 16 V"K' ^V"Z ' 
+/3,o,3 -y^^-8^ + -F +8F Z 



, V" K' V" Z V" Z'^P 128 V"'K'(t)'' 
+/3,5,o 64^^ + 64^^ + 64 ^T^+ ^ 



N N N 

128 r^^Zc/)^ 128 V" Z' (0^) '^ 
]V ^ ~3 TV 



+4,4,0 (l6 



V" K' V" Z V" Z'6^ 32 V" K'S^ 

h 16 h 16 — H — 

AT AT N Z N 



32 32 



/80 V" K' V" Z V" Z'(f)'^ 
+/3,3,o ( — + 48 + 48 



3 



A^ 



A^ 



F'"Z>22 

+32^^ + 32^^ + 32^^ + - 



/ 512 I/--^ 64 ir' 512 , 64 , / 



/ 512 1^ m 512 3 ,_128 

^'^'^ l3Ar 3Ar 3 ^ 3 



/704 7"^ 02 16 7" if' 40 V"Z 40 7" 

+ ^3 2 2 \ \ — 



N 

80 7"' ' 0^ _ 704 ^,,3 ^ 16 ^„ ^, _ 40 ^ _ 40 ^„ ^, 

3 



3 A^ 3 A^ 
K' 

aT 3 " 3 
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~ V"'K'<P^ 



3,2,1 



64 V" K' A V" Z 8 V" Z 
+ 



3 N 3 N 3 N 
16 V'K'cl)'^ 80 

y TV 



/ 112 112 



/ 128 F"' 128 



+/5,5,o (- 128 



F" K K' 
N 



V"KZ V'KZ'd)'^ V'K'ZS^ VZ^d)^ 
-128 128 —-^-128 —-^-128- 



-128 



N N N N 

V" ZZ'{(l?f 256 V'KK'cl? 256 V" K Z ct)"^ 



N 3 N 3 N 

256 V" K Z' {(t)''f 256 V" K' Z {(t)^f 256 V" Z\(t)''f 



3 N 3 N 



N 



256 



5,2,2 



208 8 V"KZ 



N 



8 y"7v"Z'c6^ 176 V" K 



N 



416 V" K' Z 



o 



3 N 3 N 

16 V" KK'^"^ 208 

"3 

176 ,0 416 



N 



8 



^ X if' + - 7" X Z + - X Z' 

N 3 3 3 ^ 

16 



2,4,0 I TV AT 3 



/ 80 V"KK' 80 

+4,1,3 -y^v 



105 



B.2 Asymptotic expressions for K and Z 



We can also calculate expressions for V,K and Z for large 0^. These were im- 
posed as boundary conditions at some suitable endpoint. A^, and denote 
constants. 

We have V{(t)^) behaving as, 
A, (02) ('^) + 

fl v/6(7V-i)^/r+^ 1 V6(l + 77) \ .,2^^) 

V24 v^v^TV 24 ^7VV^, + ^^^/5^y ^'^ > 

K{(f)'^) behaves as. 



A, (02) ^) + 2 (^v^ {l+r])(4r]' A^' + 16 ^4 

+ 11 ^fc 7?^ ^2 ^ 48 ^3 ^^3 ^ 55 ^3 ^^2 _ 52 ^3 

- 56 ^fe2 ^3 _ ^2 ^^2 _ ^^2 ^2 _ gg ^^3 ^2 

- 144 - 703 r] Ak A,^ + 336 + 40 A,,^ r] A, 

- 432 T] A J" - 80 _ Aj" + 160 ^I^T^ 

(-4^, + 2^,r/ + ^+r7^)2(5-r/)^/2^V/') + ^^/6 
yrr^ ((l2 AT + 28 7?^ + 48 77 AT - 48 AT - 32 _ 4g ^^3 

+ (6877^^ - 2477A^^ - 12 77^ AT^ + 4477^^ + 877^ 

+ l27fNA,- 16 A^ Ak" + (3 7/3 AT ^2 + 3 77 AT + 21 77'^ A,' 
+ 69 77 ^2 + 66 772 ^^2 ^ 24 ^2 + 6 772 AT yl/) + 5 77^ 
+ 15 77 A,^ + 15 772 + 5 A J) /( ^//2 ^ 
(( 277 - 4 )^, + 77^ + ^)2)) (02)(-TT^) 
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Finally Z{(f)'^) behaves as, 

A, + |^_L (64 ri^N -65rj^ + 185 r/^ - 172 rj'^ N - 716 rj^ N 

+ 652 r)^ + 604 ry^ iV - 452 rj^ - 2524 + 2348 AT - 1184 
+ 126AN)V6/(yA',N{-2 + r)y{l + r)){5-r)f/''^ 

24 (1 + ^)3/2 (_2 + ^)2^y]4; J^*^^ 
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